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FOREWORD 


The Lewis Research Center is the National Aeronautics and Space 
Administration's principal field installation for aerospace propulsion and 
power generation research and development. Therefore, a substantial part of 
the Center's activities is devoted to progress in the technology of aircraft 
propulsion. Results of this work are published as NASA reports and as 
articles in the technical journals. In addition, an occasional technical 
conference assists us in communicating more directly with others in the 
engineering field. Accordingly, this Combustion Fundamentals Research 
Conference was held to provide a forum for selected researchers from industry, 
academe, and the Government to review results of recent and current 
Lewis-sponsored research in combustion and the fluid mechanic processes in the 
gas-turbine combustor environment. This meeting also provided an opportunity 
for technical dialog among the participants. 

Over the last few years increased emphasis has been placed on fundamental 
and generic research at Lewis with less systems development efforts. This is 
especially true in combustion research, where the study of combustion 
fundamentals has grown significantly in order to better address the perceived 
long-term technical needs of the aerospace industry. The main thrusts for 
this combustion fundamentals program area as follows: 

(1) Analytical models of combustion processes - These models characterize 
the governing physical phenomena that occur during combustion and in the fluid 
dynamic processes associated with gas turbine combustors, 

(2) Model verification experiments - These experiments provide 
benchmark-quality data to use in assessing the accuracy of analytical models 
and in identifying model deficiencies, 

(3) Fundamental combustion experiments - These experiments achieve a more 
complete and basic understanding of the fundamental aerodynamic and chemical 
processes occurring in chemically reacting flows, 

(4) Advanced numeric techniques - These techniques improve computer codes 
in terms of efficiency, numerical accuracy, and display of results. 

At present, there are several research activities under each of these 
thrusts. Most of these activities were reviewed by the respective principal 
investigators. Each presentation summarizes the progress and current status 
of the project. If more detail is desired, discussion with the principal 
investigator is invited. In most cases, previous publications are listed. 
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ANALYTICAL MODELING OF OPERATING CHARACTERISTICS OF 
PREMIXING^PREVAPORIZING FUEL-AIR MIXING PASSAGES 


0. L. Anderson 

United Technologies Research Center 
East Hartford, Connecticut 06040 


A model for predicting the distribution of liquid fuel droplets and fuel 
vapor in premixing-prevaporizing fuel-air mixing passages of the direct injection 
type is described. This model consists of three computer programs: a calculation 
of the two-dimensional or axisymmetric air flow field neglecting the effects of 
fuel; a calculation of the three-dimensional fuel droplet trajectories and 
evaporation rates in a known, moving air flow; and a calculation of fuel vapor 
diffusing into a moving three-dimensional air flow with source terns dependent 
on the droplet evaporation rates. The air flow calculation can treat compressible 
swirling flows in arbitrary ducts with arbitrary distributions of temperature and 
velocity as initial conditions. The fuel droplets are treated as individual 
particle classes each satisfying Newton's law, a heat transfer, and a mass 
transfer equation. Each particle class has a number density such that summation 
over all particle classes yields the fuel flow rate. This fuel droplet model 
treats multicomponent fuels and incorporates the physics required for the treatment 
of elastic droplet collisions, droplet shattering, droplet coalescence and 
droplet wall interactions. The vapor diffusion calculation treats three-dimensional, 
gas-phase, turbulent diffusion processes with the turbulence level determined by 
the air flow calculations and the source terms determined by the droplet evaporation 
rates . 

The analysis includes two models for the autoignition of the fuel-air mixture 
based upon the rate of formation of an important intermediate chemical species 
during the pre-ignition period. This species is produced both within the vicinity 
of the fuel droplets and throughout the diffusing fuel vapor-air mixture. Since 
chemical reaction rates may depend upon the local mixture temperature, the local 
mixture temperature is adjusted for the effect of fuel evaporation. One of these 
autoignition models is based upon a global chemical reaction rate in which fuel 
decomposes to form the intermediate species, ethene. Additional chemical reactions 
result in the oxidation of ethene to carbon dioxide and water vapor. The maximum 
ethene concentration occurs at approximately the same instant in time in which 
a sudden rise in gas temperature occurs. The second model is based upon the 
global chemical reaction rate governi,ng the production of an unknown, but important, 
intermediate species; autoignition occurs where the concentration of this species 
reaches a critical value determine.d by adjusting constants in the global reaction 
rate expression to give agreement with representative, experimentally-obtained, 
autoignition times. 


1 



Model development has been completed, the computer codes have been delivered 
to NASA, and the final report has been completed. The final report includes a 
literature survey conducted to identify the important characteristics that 
influence the occurrence of autoignition and to aid in the development of auto- 
ignition models. 


OBJECTIVES 


• Develop computer codes to predict 3D distribution of 
liquid and vapor fuel in moving gas stream 

• Develop an autoignition model for fuel-air mixtures 
in LPP passages 


COMPUTATIONAL SEQUENCE 

• Calculate 2D air flow distribution 

• Calculate 3D fuel droplet trajectories and sources 
of fuel vapor and autoignition species from droplets 

• Calculate 3D diffusion of fuel vapor 

• Calculate 3D diffusion of autoignition species from 
droplet and vapor sources 

• Apply autoignition criterion 
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CALCULATION OF AIR FLOW 

(UTRC ^DD code) 

• Two dimensional or axisymmetric/swirling 

• Arbitrary duct contours 

• Arbitrary inlet flow conditions 

• Arbitrary perfect gas 

• Turbuience model 

• Mixing length free stream turbulence 

• Van Driest wall turbulence 

• Streamline coordinates-conformal mapping 

• Implicit forward marching numerical integration 


CALCULATION OF SPECIES 

(UTRC VAPDIF code) 


Dt 


~ 2 

V C = Sg + Sj 


• 3D binary diffusion equa'ion 

• Fuel vapor into air 

• Autoignition species into air 


• Source terms 

• Fuel vapor from droplet evaporation 

• Fuel vapor from evaporation at a wall 

• Autoignition species from droplet evaporation 

• Autoignition species from diffusing fuel vapor 

• Turbulent eddy diffusivity 
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CALCULATION OF FUEL 
DROPLET EVAPORATION 

(UTRC PTRAK code) 

• Air flow field known 

• Droplets treated as particle classes 

• (3) momentum equations 

• Mass transfer equation 

• Heat transfer equation 

• Fuel liquid and vapor thermodynamics 

• Multicomponent fuel blends 

• Droplet shattering, coalescence, and 
well interactions 


AUTOIGNITION MODEL I 

Hautman, Dryer, Shug and Glaaaman 
• Fuel (CnH2n+2) 

K Y°1 

— A I- A — 


dt 




Ethene (C2H2) 
dXc 


dt 



Time 


^2 y "^2 

O 2 F , E 
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AUTOiGNITION MODEL 1 - (Cont.) 

• Provides general analytical frame>work 

• Requires strong interaction between fuel vapor 
and autoignition species 

• Requires calibration at LPP conditions 


AUTOIGNITION MODEL II 


• Autoignition species (?) 
dlXe/Xg*) 


dt 

Criterion 


= KXfXo2<#>^tV 


f'^ dCXg/XE*) 

L dt = 1 

o dt 


• Correlation 

X = T(P,T,<#>,XfXo^...| 

• Simple to apply 

• Difficult to calibrate 
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NUMERICAL SIMULATION FOR DROPLET COMBUSTION USING LAGRANGIAN HYDRODYNAMICS 


Martin J. Fritts, David E. Pyfe, and ELaine S. Oran 
Laboratory for Computational Physics 
Ifeval Research laboratory 
Washington, DC 20375 

A predictive model of spray combustion nust incorporate models for the wide 
variety of physical environments in a practical combustor. In regions where 
droplets are closely spaced, combustion resembles a diffusion flame-, where they 
are well-separated, an envelope or wake flame results. The relative velocity 
field between the fuel droplets and oxidizer influences boundary layer develop- 
ment about the droplet, recirculating flow patterns, and droplet shape and 
stability. A model mist encorpass these interacting tenporal and spatial 
effects as well as corcplicated combustor boundaries. The objective of the 
current work is to develop the triangular gridding method for describing the 
individual and collective properties of vaporizing and burning fuel droplets. 

Our approach to this problem has been to modify the basic two dimensional 
Lagrangian model [l] to simulate flows in and about fuel droplets [2]. This 
general-connectivity triangular grid permits accurate representation of 
boundaries as well as material surfaces and interfaces. It also allows variable 
resolution through the insertion of new cells as required to maintain accuracy. 
The finite-difference operators for divergence, curl and gradient are 
constructed to exactly reflect the properties of the continuum operators. The 
construction assures conservation of vorticity and mass and provides a 
determination of the local grid connectivity based on convergence criteria for 
the solution of Poisson’s Equation. Extensions of the model have been made in 
the same spirit: finite-difference operators conform to the continuum limit and 

physical quantities are conserved properly. An important factor to note is that 
development of this method is entirely original work; there is no basic lore to 
fall back on when something goes wrong, and algorithms for a particular type of 
term mist be devised and the best one chosen. 

To date the basic hydrodynamic code incorporates algorithms which allow us 
to include the effects of surface tension and viscosity. These algorithms have 
been incorporated and tested extensively as described below. We have also 
devised algorithms for including the effects of compressibility for subsonic 
flow and for incorporating the effects of thermal conductivity. While these 
latter effects are being tested and incorporated, we will be developing the 
algorithms for vaporization and molecular diffusion. The final steps will allow 
us to describe a burning droplet. 

The test problems performed have included sirailations of inconpressible 
flows about droplets in which the density ratios of droplet to background mate- 
rial have been 2:1, 10:1, and most recently 800:1. The two lower density ratios 
were used to test the various algorithms. The latest tests are aimed at model- 
ling kerosene in air and have incoporated the egressions for surface tension 
and viscosity. 

As an exanple of the way in which development of algorithms has proceeded, 
we describe the latest development in an algorithm for viscosity. Our viscosity 
algorithm originally ewressed the change in the vorticity, at a grid point 
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due to viscosity by d^/dt = where v is the viscosity coefficient. All 

triangle velocities about a vertex contributed equally to the change in vorti- 
city at the central grid point. Although the algorithm produced the correct 
spreading rates for the test case of a shear profile, we found that it only did 
so for fairly regular grid geometries because of the ambiguity in determining 
how the changes in vorticity are translated to velocity changes for different 
grid geometries. For an arbitrary grid, a more complete prescription was neces- 
sary and a method was developed in which V^V is a triangle-centered 
quantity. 

This new triangle -centered algorithm was tested in a calculation of the 
spreading of a shear layer of initially zero thickness. The way in which the 
velocity distribution across this layer evolves and the growth of the width of 
the layer are known quantities which may be coiripared to the results of calcula- 
tions. We found that the calculated layer width agreed exactly with the theory 
and the shear layer velocities were correct over the entire mesh. The 00115 ) 0 - 
nents of velocity perpendicular to the shear layer remained zero indicating that 
the algorithm worked well even for the distorted grid used in the test prob- 
lems. 


In order to test our algorithm for surface tension, we performed calcula- 
tions of droplets which oscillate under the effects of surface tension. The 
results of sinulations could be conpared to the linear theory for small anpli- 
tude oscillations on cylindrical Jets, as first given by Rayleigh. We extended 
this theory to predict droplet oscillation frequencies for a droplet in a back- 
ground gas of finite density. 

In the numerical calculations we studied an n = 2 oscillation for a droplet 
density of 2 g/cm^ and an external fluid density of 1 g/cm^. From the calcula- 
tions shown in Figure 1 , we find that the n\americal oscillation period is 
approximately 1.25xl0“^s whereas the theoretical period is 1.13><10"^s. Most of 
the small discrepancy between the numerical and theoretical results can be 
explained by the finite grid s^'^^cing. However, given Rayleigh’s experience with 
large anplitude oscillations, it is reasonable to expect our conputational 
period to differ somewhat from that given by the linear theory. 

Figures 2 and 3 show other early test calculations. Figure 2 illustrates 
the case for which the density ratio is 2:1 and there is no surface tension or 
viscosity present. We see that a recirculation zone forms early in the calcula- 
tion, compressing the droplet in the direction parallel to the flow. Flow 
within the droplet is initiated by this conpression in a direction normal to the 
external flow. The bulges formed at the top and bottom of the distorted droplet 
are pulled around the recirculation zone by the shear flow which is at a maximum 
at these points. The internal droplet flow is therefore driven by the corrpres- 
sion set up between the front and rear stagnation points and by the high shear 
flow which extends around the top and bottom of the droplet and recirculation 
zone. Ihe interaction of the droplet back onto the external flow occurs pri- 
marily throiigh the enlarged cross-sectional area of the droplet which increases 
the size of the recirculation zone. Eventually the droplet is squeezed into a 
thin layer coating the recirculating zone. The thinned film then shatters into 
several smaller pieces, first at the rear of the droplet and later in the more 
laminar flow toward the front on the droplet. 


8 



Figure 3 shows the results of a calculation with surface tension for the 
same initial conditions as used in the calculation without surface tension 
(Figure 2). As in the case without surface tension, the internal droplet flow 
is driven by conpression parallel to the external flow and is initially norml 
to the external flow. A recirculation zone is formed in the wake of the 
conpressed droplet and the droplet is deformed as it is swept outward and back- 
ward by the external flow both outside and inside the recirculation zone. 
However, the presence of surface tension provides counteracting forces at 
regions of high curvat\ire. Such forces at the sides and rear of the droplet are 
sufficient to stop the droplet from thinning around the recirculation zone. 

The presentation will summarize the calculations shown above and then 
proceed to describe the more recent calculations at an 800;1 density ratio. 
Introducing this large ratio initially caused several problems in defining 
velocities when cells near the droplet boundary were divided. Ihis problem was 
due to finite resolution effects at the boundary and has been fixed. Current 
simulations with viscosity should give us a realistic picture of flows in the 
droplet itself. 

Finally, we describe the new algorithms to be tested for subsonic conpres- 
sibility and thermal conductivity. A new parametric representation of taut 
splines will also be presented. Ihis extension was necessary to achieve a 
”wiggle-free” spline fit in smooth regions of a droplet interface which were 
near surface discontinuties. 
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STRUCTURE OF EVAPORATING AND COMBUSTING SPRAYS: 
MEASUREMENTS AND PREDICTIONS* 


J-S. Shuen, A.S.P. Solomon and G. M. Faeth 
The Department of Mechanical Engineering 
The Pennsylvania State University 
University Park, Pennsylvania 16802 


Summary 

Results obtained during the first year of an investigation of the 
structure of sprays are briefly described. Further details may be found 
in Refs. 1-3. 

The investigation involves both experimentation and analysis. Experi- 
mental objectives are to complete measurements of the structure of nonevapo- 
rating, evaporating and combusting .sprays for sufficiently well-defined 
boundary conditions to allow evaluation of models of these processes (Fig. 1), 
Analytical objectives are to begin model evaluation using both existing and 
the new data (Fig. 2). The results of the investigation have application to 
the development of rational design methods for aircraft combustion chambers 
and other devices involving spray combustion (Fig. 3). 

Major assumptions for the models are summarized in Fig. 4. The con- 
tinuous phase is treated using a k-e-g model of turbulence originally pro- 
posed by Lockwood and Naguib [4], which has been extensively calibrated for 
noncombusting and combusting single-phase flows during earlier work in this 
laboratory [ 5-7 ] . 

Three methods for treating the discrete phase are being considered: 

(1) a locally homogeneous flow (LHF) model, (2) a deterministic separated 
flow (DSF) model, and (3) a stochastic separated flow (SSF) model. The 
main properties of these models are summarized in Figs. 5-7. Infinitely 
fast interphase transport rates and local thermodynamic equilibrium are 
assumed for the LHF model (Fig. 5) — implying that both phases have the same 
temperature and velocity at each point in the flow. LHF models provide a 
useful limit for infinitely small particles or drops, but generally over- 
estimate the rate of development of practical sprays [5-7]. DSF models 
(Fig. 6) provide for finite interphase transport rates, but assume that 
interphase transport can be found by ignoring effects of turbulent fluctu- 
ations. Most spray models reported to date employ this approximation. The 
present SSF model (Fig. 7) adapts an approach originally proposed by Gosman 
and loannides [8]. In this case, particles or drops are assumed to interact 
with a succession of, turbulent eddies whose properties are determined by 
random sampling — given mean and fluctuating properties of the flow from the 
k-c-g model calculations. This involves computation of a statistically 
significant number of particle trajectories using Monte Carlo techniques. 


* 


NASA Grant No. NAG 3-190 with R.' Tacina of Lewis Research Center as 
NASA Scientific Officer. 
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Initial model evaluation employed data for dilute particle-laden jets — 
to avoid complications due to particle coalescence, particle collisions and 
polydisperse particle flows. Only a sample of the results is given here, 
cf. Refs. [1-3] for complete findings. 

The prescription for eddy properties used in the SSF model was cali- 
brated using theoretical results of Hinze [9] — similar to Gosman and 
loannides [8]. This analysis was for the dispersion of inf initely-small 
particles in a homogeneous and isotropic turbulent flow. The comparison 
between SSF predictions and the analytical result is illustrated in Fig, 8, 

Good agreement was achieved — fixing methods for estimating particle inter- 
actions with eddies. The SSF model predictions are compared with measure- 
ments of particle dispersion in a duct flow, reported by Snyder and Lumley 
[10], in Fig. 9 , The model is seen to provide encouraging predictions of 
effects of particle properties on rates of particle dispersion. 

The remaining comparisons between predictions and measurements consider 
particle-laden jets. LHF, DSF and SSF model predictions are compared with 
the measurements of Yuu et al, [11] in Fig. 10. The LHF and DSF models 
over- and under-estimate particle dispersion, while the SSF model is in good 
agreement with measurements. Furthermore, the DSF model indicates that the 
particles tend to concentrate near the centerline as axial distance increases — 
which is not observed. Comparison of the models with measurements of McComb 
and Salih [12,13], Laats and Frishman [14,15] and Levy and Lockwood [16] 
continues in Figs. 11-13. The SSF model yields satisfactory predictions for 
these flows, aside from possible effects of turbulence modulation and turbu- 
lence generation at high particle mass loadings [2,3]. This evaluation is 
not adequately definitive, however, due to uncertainties in initial conditions 
for the existing particle-laden jet data [1-3], 

Evaluation of the models is continuing using data from the present in- 
vestigation. A sketch of the test apparatus being used for noncombusting 
sprays is illustrated in Fig. 14. An air-atomizing injector sprays vertically 
downward along the centerline of a traversible screened enclosure. Experi- 
mental methods are summarized in Fig. 15. All techniques have been used in 
work to date, aside from. the LDA-visibility method for drop size and velocity 
measurements — which requires a different optical geometry and is being de- 
ferred until other measurements are complete. 

Experimental methods were established by satisfactory measurements of 
the properties of air jets, formed by the injector, with earlier work [5,6]. 
Tests were then conducted in two nonevaporating sprays having SMD of 87 and 
30 ym. LHF and SSF model predictions are compared with measurements of mean 
gas velocity and mean liquid flux, along the spray axis, in Figs. 16 and 17. 
There is no fundamental limitation in the use of the LHF model in dense 
regions of the spray; therefore, these predictions extend from the injector 
exit. The separated flow models, however, are limited to dilute regions of 
the spray; therefore, these calculations begin at x/d = 50 — where adequate 
initial conditions were available from the measurements. The prediction of 
the LHF model improves for the more finely atomized spray. Case 1, but is 
not very satisfactory. In contrast, the SSF model provides reasonably good 
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predictions of these measurements. Predicted and measured radial profiles 
of liquid mass flux are plotted in Fig. 18. Similar to the results for 
particle-laden jets, the DSF model yields excessive concentrations of the 
dispersed phase near the axis, since turbulent particle diffusion is ignored. 
The SSF model provides fair predictions of particle spread. However, in 
this case, the LHF model underestimates spread rates I This effect is well- 
known, involving enhanced spreading of particle-laden flows by turbulent 
diffusion for a certain range of particle inertial properties. The fact 
that the SSF model correctly predicts this trend is very encouraging. Pre- 
dicted and measured turbulence kinetic energy are illustrated in Fig. 19. 
While the LHF model underestimates the magnitude of k and the width of the 
flow, the SSF model provides satisfactory predictions. This, too, is 
encouraging, since predictions of k are an important element in estimating 
eddy properties for the SSF model. Good predictions of Reynolds stress 
were also obtained with the SSF model. Reynold stress depends on z pre- 
dictions — suggesting that this aspect of the eddy prescription is also 
adequate. 

The conclusions, to date, are summarized in Fig. 20. The main con- 
clusion is that the SSF model provides encouraging predictions for these 
multiphase flows — with minimal added empiricism. It will be most interesting 
to examine this methodology for evaporating and combusting sprays — where 
effects of concentration fluctuations must be considered along with velocity 
fluctuations. The main limitation of the evaluation, thus far, is adequate 
specification of initial conditions for the present data base. Tests during 
the next report period are designed to eliminate this deficiency (Fig. 21). 
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pp. 333-339. 
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Objectives:’ Complete measurements of spray structure suitable 
for evaluation of models. 

Configuration: Axisymmetric spray or particle-laden jet in a 

quiescent environment. 

Specific Cases: i. Air jet (calibration). 

II. Particle-laden jet (sand particles), 
hi. Non-evaporating spray (low-volatility oil). 
IV. Evaporating spray (Freon-11). 

V. Combusting spray (h-Pentane). 


Figure 1.- Experimental Objectives. 


Objectives; Complete evaluation of typical models using both 
existing and new data. 

Models: i. Locally homogeneous flow (LHF)— Infinitely 

fast interphase transport rates. 

ii. Deterministic separated flow (DSF)— Finite 
interphase transport rates considering 
particle response to mean motion. 

III. Stochastic separated flow (SSF)— Finite 
interphase transport rates but particles 
respond to individual eddies. 

Evaluation: Gas jets^ solid-particle-laden jets, nonevaporating 
sprays, evaporating sprays and combusting sprays. 


Figure 2.- Analytical Objectives. 
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Potential contributions of the results are: 

1. The development of rational design methods for 

AIRCRAFT combustion CHAMBERS^ GAS-TURBINE 
COMBUSTORSy FURNACESy STRAT I F I ED-CHARGE I.C. 
ENGINES AND DiESEL ENGINES, 

2. Improved understanding of the properties of 

TURBULENT PART I CLE/DROP-LADEN FLOWS, 


Figure 3.- Applications of the investigation. 


1. Steady^ AXisvMMETRiCy large Reynolds NUMBERy boundary-layer flow, 

2. k-c-g TURBULENCE MODEL WHICH IS WELL-CALIBRATED FOR NONCOMBUSTING 
AND COMBUSTING SINGLE-PHASE JETS, 

3. Negligible kinetic energy and viscous dissipation of mean flow 
AND radiation. 

Exchange coefficients of all species and heat identical. 

5, Continuous-phase is in local thermodynamic equilibrium. 

6. DSF AND SSF ONLY; Dilute -particulate flow so effects of turbulence 

GENERATION AND DISSIPATION BY PARTICLESy PARTICLE COLLISIONSy AND 
ADJACENT-PARTICLE DISTURBANCES OF INTERPHASE TRANSPORT RATES ARE 
NEGLIGIBLE. 


Figure 4.- Major assumptions of the models. 
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Assumption: Infinitely fast interphase transport rates, i,e., 

PARTICLE AND CONTINUOUS PHASE VELOCITIES AND 
temperatures are IDENTICAL AND LOCAL THERMODYNAMIC 
EQUILIBRIUM INCLUDES BOTH PHASES. 



Notes: i. 

II. 

in. 


Correct limit for infinitely small particles. 

Maximum particle response to turbulent fluctuations. 
Computation equivalent to single-phase flow. 


Figure 5.- Properties of the LHF model. 


Assumption: Finite interphase transport rates with particles 
responding to mean motion. 



pi 


pi 


^ U, 


ip 


T 

^G. 


Notes: i. Particle dispersion ignored— only valid for 

"large" particles. 

II. Typical approach in current spray models. 

III. Eulerian calculation for continuous phase with 

DISTRIBUTED SOURCE TERMS FROM PARTICLE INTERACTIONS. 

IV. LaGRANGIAN calculation of PARTICLE TRAJECTORIES. 

Figure 6.- Properties of the DSF model. 
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Assumption: Finite interphase transport rates with particles 

INTERACTING WITH A SUCCESSION OF INDIVIDUAL EDDIES 
WHOSE PROPERTIES ARE FOUND BY RANDOM SAMPLING OF 
LOCAL TURBULENCE PROPERTIES. 



, (2k/3)’''^) 

t|pUF(?, g’ 
= G.jpDF (7 


Notes: i. Maximum particle displacement and time of eddy 

• p 3/4 i3/2y . _ I 

INTERACTION ARE: *-e ' ^ ' Lg/(2k/3) 

II. Provides predictions of fluctuating particle 

PROPERTIES AND TURBULENT PARTICLE DISPERSION. 


III. Computation similar to DSF model--Monte Carlo 

TECHNIQUE TO FIND PARTICLE TRAJECTORIES. 

Figure 7.- Properties of the SSF model. 



Figure 8.- SSF predictions for the dis- 
persion of infinitely small particles 
in a homogeneous isotropic flow (ana- 
lytical results from Hinze [9]). 
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0 100 200 300 400 

t(ms) 


Figure 9.-. Predicted and measured particle dispersion in a 
uniform grid-generated turbulent flow (measurements from 
Snyder and Lumley [lO]). 



Figure 10.- Comparison of LHF, DSF, and SSF 
predictions of particle dispersion with 
the measurements of Yuu, et al. [ll]. 
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r/x 

Figure 13.- Predicted and measured mean 
and fluctuating particle velocities 
(data from Levy and Lockwood [16]). 



Figure 14.- Sketch of the experimental 
apparatus. 


23 



Measurement 


Technique 


Mean and fluctuating gas and 

PARTICLE VELOCITIES 

Drop (particle)— size distributions 

Drop-size distributions 

Drop size and velocity 
Mean liquid (particle) flux 

Mean temperature 
Mean composition 


Laser-Doppler anemometer (LDA) 

Fraunhofer diffraction (monitor- 
ing nonevaporating flows only) 

Slide impaction 

LDA-visibility method 

Isokinetic sampling and filtering 
(nonevaporating flows only) 

Shielded, fine-wire thermocouple 

Isokinetic sampling and analysis 
with gas chromatograph 


Figure 15.- Summary of experimental methods. 



x/d 


Figure 16.- Predicted and measured mean gas 
velocities along the axis of nonevaporating 
sprays . 
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Figure 17. Predicted and measured mean liquid 
mass flux along the axis of nonevaporating 
sprays . 











1. Model Performance: 

I. LHF AND DSF MODELS ARE VALID AT LIMITS OF SMALL AND LARGE 
PARTICLES^ BUT DATA BASE (AND PROBABLY SPRAYS IN GENERAL) 
INCLUDED FEW RESULTS AT THESE CONDITIONS. 

II. SSF MODEL (PARTICULARLY THE NSSF VERSION) YIELDED ENCOURAGING 
RESULTS— WITH MINIMAL EMPIRICISM. 

2. Decisive model evaluation was not achieved due to inadequate 

SPECIFICATION OF INITIAL CONDITIONS FOR MOST OF THE DATA BASE. 

3. Effects attributable to turbulence generation and modulation 

WERE OBSERVED^ HOWEVER^ THE EFFECTS WERE SMALL SINCE. DATA BASE 
WAS LIMITED TO DILUTE SPRAYS (VOID FRACTION > 99%). 


Figure 20.- Conclusions of the investigation* 


Task 

Status 

Model Development 

Completed for noncombusting 
flows. Combusting flows next. 

Measurements 


Single-Phase Jets 

Completed. 

Particle-Laden Jets 

Completed, aside from mean and 
fluctuating particle velocities. 

Nonevaporating Sprays 

Completed^ aside from drop size 
and velocity measurements. 

Evaporating Sprays 

Completed^ aside from mean and 
fluctuating gas velocities (in 
progress) and drop size and 
VELOCITY measurements. 

Combusting Sprays 

Next. 


Figure 21.- Status of the investigation. 
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FUEL SPRAY DIAGNOSTICS 


Francis M. Humenik and >fa.ria A. Bosque 
NASA Lewis Research Center 
Cleveland, Ohio 44135 


BACKGROUND : 

The two most significant parameters which characterize fuel 
spray combustion are the droplet size and the size 
distribution. Since the basic interactions between fluid 
mechanics and the combustion process are not well understood, a 
great deal of research, both experimental and theoretical, has 
addressed itself to the problem of spray combustion. Turbulent 
fuel and air motions influence the chemical reactions by 
increasing the oxygen supply to the fuel. The relative 
velocities between the fuel's gas phase and liquid droplets 
affect the evaporation, burning rate and the pollutant 
formation. Thus, it is necessary to measure the size and 
velocity distributions of the spray droplets and the evolution 
of these distributions with the flow. Such measurements will 
enable researchers to get a better picture of spray combustion 
and in turn help in both fuel nozzle research and developing 
theoretical models of spray combustion. 

The Droplet Sizing Interferometer has promised to be very 
successful in producing droplet size and velocity and 
size-velocity correlations measurements under a variety of 
spray conditions. By using off-axis large angle light scatter 
detection, the measurement region has been substantially 
reduced and spatial resolution was significantly improved. 

Using the light scattered by reflection and refraction made the 
measurement of droplet size and velocity possible in dense 
sprays. Such measurements are important to enable comparisons 
of different fuel injectors and to characterize the combustion 
processes in gas turbine engines. 


B. RESEARCH OBJECTIVES: 


Using water as a fuel substitute, severa:! fuel nozzles will 
be selected for measurement of droplet size and velocity at 
representative test conditions. In addition, the turbulent 
interactions between flowing air streams and simulated fuel 
sprays will be studied. Later experiments will investigate 
simple flow fields using small particles for air flow seeding 
to track the flow characteristics. 

C. TECHNICAL SIGNIFICANCE: 


These experiments will provide the fundamental experimental 
data base for turbulent flow mixing models and lead to better 
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prediction of more complex turbulent chemical reacting flows. 
The experiments will also be useful for analytical application 
to combustor design and provide a better fundamental under- 
standing of the combustion process. 

D. APPROACH : 

Several experimental configurations are currently being 
considered for laser diagnostic investigation. The initial 
experimental configuration will consist of a simple mounting 
device for fuel nozzle characterization tests. The device will 
have manual traversing capability to get several probe volume 
positions for mapping the fuel nozzle flow distribution and 
variations with test conditions. Another experimental 
configuration would utilize a fuel nozzle/swirler combination 
to study simulated fuel/air mixing characteristics. A possible 
simple turbulent mixing study would utilize a flexible arrange- 
ment of clear tubular lucite sections in a configuration to 
permit flow seeding and adjustable positioning features to 
permit manually scanning of the turbulent flow field. Details 
of these various configurations will be developed, pro- 
gressing with additional experience and developing expertise 
with the laser diagnostic equipment. 

V/ith each experimental configuration, the Droplet Sizing 
Interferometer apparatus will be used to obtain a series of 
measurements to fully map droplet size and velocity distri- 
butions. Two receiver units will be positioned off-axis at a 
collection angle selected to cover the expected size range with 
acceptable signal/noise ratio. Each receiver unit contains the 
collection optics, photomultiplier tube assembly, and high 
voltage power supply. A visibility and Doppler signal 
processor accepts signals from the PMT and simultaneously 
processes each signal for visibility and Doppler period which 
correspond to size and velocity. Data from the visibility 
processor are accepted and displayed by the Data Management 
System (DMS). 


E. DATA SYSTEM DESCRIPTION: 


The Droplet Sizing Interferometer (DSI) consists of 13 
major components which include a transmitter unit, two receiver 
units, two signal processors, two data management systems, two 
Bragg cell systems, two printer/plotters, a laser, power 
supply, and a color monitor. The transmitter unit includes a 
0.5 watt Argon-Ion Laser, beam steering devices, frequency 
shifting capability, beam splitting and focusing components. 

The system is a two-color, two component system. Two 
independent, orthogonal measurements of size and velocity 
components can be made simultaneously. 
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The Data Management System (DMS) includes a 6502-based 
micro-computer, dual disk drives, video display, and I/O 
hardware to interface the visibility processor. A printer/- 
plotter is included to provide hard copy records. The DMS is 
designed to accept data from the visibility processor and 
display the data in real-time in histogram form on the video 
monitor. The data can be stored and retrieved from mini-floppy 
disks. The DMS also handles data reduction. The software 
package provides a variety of mean diameter calculations, 
velocity calculations, and data plots. 


F. STATUS: 

Initial training and practice sessions have been conducted 
during the past year with simple laboratory spray experiments. 

The test cell for formal spray diagnostic experiments is 
now operational; experimental spray configurations are in the 
process of being set-up. The test facility is presently 
equipped with 2 fixed optical tables, experiment breadboards, 
and other optical accessories. A traversing optical table 
system is on order with delivery expected in 4-6 months. 

Two contracts have been awarded for further development and 
improvements to the system. These contracts promise to extend 
the size range capability, reduce beam alignment difficulties, 
and reduce the system sensitivity to laser beam quality and 
differences in the relative intensity of the beams. 

Data obtained to date have been limited to verification 
experiments with a monodisperse droplet generator and system 
familiarity experiments. 


G. REFERENCES; 


Bachalo, William D. : Method for Measuring the Size and 

Velocity of Spheres by Dual-Beam Light -Scatter Interferometry; 
Applied Optics, Vol. 19, P. 363, Feb. 1980. 
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LASER DIAGNOSnC PROGRAM 


OBJECTIVE: 


SIGNIFICANCE: 


PURPOSE • 


METHOD 


YIELD ■ 


• MEASURE SPRAY DROPLET SIZES AND VELOCITIES 
SIMULTANEOUSLY 

• MEASURE TWO COMPONENT TURBULENT FLOW PARAMETERS 

• CHARACTERIZE FLOW MIXING PROCESSES 

• EVALUATE FUEL INJECTOR FLOW DISTRIBUTION PATTERNS 

• FUNDAMENTAL DATA BASE FOR TURBULENCE MODEL DEVELOPMENT 

• EVALUATE TURBULENCE INTENSITY AND REACTION CONDITIONS 
IN COMBUSTION ZONES 


LASER DIAGNOSTIC FACILITY 


CONDUCT FUNDAMENTAL COMBUSTION PROCESS RESEARCH USING 
DIAGNOSTIC MEASUREMENTS FROM A UNIQUE LASER INTER- 
FEROMETER APPARATUS. 


NON-INTRUSIVE MEASUREMENT PROBE FORMED BY COHERENT 
LASER BEAMS THAT PRODUCE INTERFERENCE FRINGE BANDS. 

DROPLETS/PARTICLES GENERATE LIGHT SCATTER PATTERNS 
WHICH ARE DETECTED WITH PMT. 


FAST REAL TIME DATA ANALYSIS OF FLOW CHARACTERISTICS. 
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INTERFEROMETER APPARATUS 


DESCRIPTION • ARGON- ION LASER 
LEXEL MODEL 85-. 5 
0.5 WATT POWER 

• TRANSMITTER TWO RECEIVER SYSTEM 
SPFCTRON LABS MODEL AOOO 
TWO-COLOR. TWO-COMPONENTS WITH ORTHOGONAL 
PAIRS OF BEAMS 


CAPABILITIES • DROPLET SIZING 

OFF-AXIS DETECTION 
3 TO 3000 MICROMETERS 
THREE DECADE RANGES 

• TWO VELOCITY COMPONENTS 
UP TO TOO M/SEC 
ELEVEN RANGES 
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LASER DIAGNOSTICS FACILITY 


BLUE 


TRANSMITTER 


Y HORIZONTAL 
VELOCITY 
COMPONENT 


LASER 


PROBE VOLUME 


LEXEL HVPS 


BRAQQ CELL 
DRIVER 


VERTICAL 

VELOCITY 

COMPONENT 


VISIBILITY 

PROCESSOR 


VISIBILITY 

PROCESSOR 


COLOR 

MONITOR 


DATA 

MANAGEMENT 

SYSTEM 


(MASTER) 


PRINTER 


keyboard 


KEYBOARD 


DATA 

MANAGEMENT 

SYSTEM 



(SLAVE) 


PRINTER 











LASER DIAGNOSTICS PROGRAM SCHEDULE 


.jITIMli- 








Development of two measurement techniques 


-The Visibillty/Intensity Technique (V/l) 

-The Dual Beam Maximum intensity Thechnique (IMAX) 


Limitations of the Droplet Sizing Interferometer (DSI) 
in practical enviroments : 


-Particles prior to crossover alter frln*jo pattern 

-Multiple particles in probe volume. . 

-Visibility Is a function of where the particle goes 
through the probe volume. 


PARTICLE SIZING INTERFEROMETRY 



Size Histogram of a Monodisporsed 
Spray in an Undisturbed Probe 
Volume (Visibility Technique) 


MONODISPERSE 



Size Histogram of a Hoiu^d ispersed 
Spray (Visibility lecimiqiie) . 
Laser Blocked by Spray. 
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THE VISIBILITY/INTENSITY TECHNIQUE (V/l) 


-Makes use of the absolute Intensity of the scattered light to provide 
a criterion for signal validation. 

-Method to establish limits for every measured visibility: 

droplets that produce must have a must scatter light 

certain visibility given size with a given Intensity 



Visibility and Intensity 
of a droplet traveling 
thru the middle of a 
probe volume. 


THE DUAL BEAM 

MAXIMUM INTENSITY TECHNIQUE (IMAX) 

-Eliminates the ambiguity of the Gaussian beam Intensity distribution 
and provides a direct relationship between scattered light Intensity 
and droplet size. 

-Method: 

A small beam crosses through the middle of a big beam to produce 
fringes in a region of uniform intensity In the big beam. 

OPTICAL configuration: 
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POLLUTANT FORMATION IN MONODISPERSE 
FUEL SPRAY COMBUSTION 


N. P. Cernansky and H. Sarv 
Mechanical Engineering Department 
Drexel University 
Philadelphia, Pennsylvania 19104 

The combustion of liquid sprays represents an extremely important 
class of combustion processes. In the transition region, encompassing 
droplet sizes in the range of 25-80 ym diameter, the mixing and evap- 
oration processes are both Incomplete at the flame front and burning 
occurs in a combined diffusive and premixed fashion. Under these condi- 
tions, the relative Importance of heterogeneous and homogeneous effects 
in dominating the combustion process is switched and gives rise to a 
number of interesting phenomena. For example, maxima in burning velo- 
city, extended flammability limits, minima in ignition energy, and 
minima in NO^ emissions have all been observed and reported. However, 
the actual mechanism and Important physical processes controlling the 
anomalous NO behavior in such systems are not clear at this time and 
need to be explored. 

Consequently, in November 1979, a study of NO^ formation in mono- 
disperse spray combustion was initiated with the following specific 
objectives: 

1) To quantitatively determine the effect of droplet size, number 
density, etc. on NO^ formation in monodisperse fuel spray com- 
bustion; and 

2) To isolate the Important physical and chemical phenomena in NO^ 
formation in these combustion systems. 

The experimental facility developed for this study has been improved 
at different stages. In its current configuration the set up consists of 
a Berglund-Liu Vibrating Orifice Monodisperse Aerosol Generator which 
produces monosized droplets within 1% of the mean droplet diameter. The 
spray of droplets Is dispersed, diluted and burned in a one dimensional 
flame) stabilized on a water cooled screen flsune holder. The flame is 
surrounded by a long pyrex tube to avoid outside air entrainment. A 
combination gas sampling/thermocouple probe is used for NO^ concentration 
and temperature measurements. 
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Measurements have been made over the monodisperse operating range of 
the system encompassing droplet diameters from 36 to 70 ym and equivalence 
ratios from 0.8 to 1.2# Radial profiles confirm the one dimensionality of 
the combustion system. Prevaporized and premixed conditions were examined, 
as well, to give the small droplet size limit. Different hydrocarbon fuels 
such as isopropanol, methanol, n-heptane, n-octane, and isox>ctane have been 
used for detailed experimental measurements. 

The experimental results indicate that both NO and NO decrease with 
decreasing droplet diameter in the spray, reaching a minimum around 48— 

58 ym for the fuels tested, and then increase again to a fairly constant 
value with further decrease in droplet diameter. The apparent behavior 
is due to increased droplet interactions as droplet size is reduced, 
causing a local oxygen depletion around burning droplets which results 
in a subsequent reduction in flame temperature and NO . Further reduc- 
tion in droplet diameter increases NO^, ultimately reaching the constant 
premixed value. 

A shift in NO minima was observed for the fuels tested; these shifts 
X 

were consistent with differences in the computed evaporation parameters. 
Calculations indicated that, relative to .isopropanol, fuels with a higher 
evaporation ratio such as methanol require larger droplets initially to 
produce the same level of prevaporation. On the other hand, n-octane and 
n-heptane, with a lower evaporation ratio, achieve the same degree of pre- 
vaporization at smaller initial droplet sizes. 

The extent of prevaporization and its effects on NO^ formation was 
further examined by preheating the air, resulting in enhanced evaporation 
of droplets. As expected, with more fuel in the vapor form, the minimum 
NO point was shifted towards larger initial droplet sizes. 

X 

Multicomponent fuel studies and synthetic oxidizer experiments are 
currently underway in order to further elucidate and quantify the impor- 
tant processes. 
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Background 


• Spray Burn i ng W i despread 

- 30-50% of Total Energy Consumption 


• Modes of Spray Combustion 

- Premixed 

- Diffusive 

- Trans i t ion 


•Transition Region Effects 

- Increased Flame Speeds 

- Broadened Flammability Limits 

- Lower Ignition Requirements 

- Reduced Emissions 


Previous Work Demonstrated an 


Anomalous NOx Behavior 


• Droplet Size Affects NO^ 

• Minimum NOx Formation at an Optimum Droplet 
Diameter 

• The Optimum Droplet Diameter Depends on Fuel 
Properties and Evaporation Character i st i cs 
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To Gas Analyzers 

ll/rj};*— Pyrcx Chimney (26.7cm X 
Screen 3.2cm ID) 
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Upper Housing 
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AXIAL TEMPERATURE. AND NO^ PROFILES ALONG THE BURNER CENTERLINE 
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«PP-> NO, NO (pp.) 





DROPLET SIZE EFFECT USING ISOPROPANOL AT 8.382 CC/MIN 





DROPLET DIAMETER fymi 


DROPLET SIZE EFFECT USING N-HEPTANE AT 0.382 CC/MIN 
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TEMPERATURE C*C) TEMPERATURE C C> TEMPERATURE C*C> 





Computed Fuel Evaporotion Parameiers 


Fuel 

Mass Transfer 
Number ^ B 

Evaporat ion 
Constant 
Kxl0 ^ (cm^/sec) 

Evaporation 
Mass Ratio 

Fu«l Isopropanol 

Isopropanol 

0.197 

0.445 

1.000 

N-Oclane 

0.017 

0.048 

0.079 

N-Heptane 

0.089 

0.243 

0.787 

Methanol 

0.482 

0.963 

3.608 


Relative 

Droplet Evaporat 

ion Effects 

Fuel 

Initial 

Diameter, jjm 

Final Relative 
Diameter, >jm 

Isopropanol 

50 

49.0 

N-Octane 

50 

49.9 

N-Heptane 

50 

49.6 

Methano 1 

50 

46.3 
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Summary of Results and Conclusions 


e Significant Effects of Droplet Diameter on 
NOx Formation, With a Minimum NOx Operating 
Condition Independent of Aerodynamic and 
Mixing Effects 

o This Anomalous NOx Behavior is Due to The 
Transition From Diffusive to Premixed type 
of Burning 

o Shifts in NOx Minimum Point ere Consistent 
With The Physical Properties and Evaporation 
Characteristics of The Fuels 
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FUEL VAPORIZATION EFFECTS 


I^ria A. Bosque 
NASA Lewis Research Center 
Qeveland, Ohio 44135 


Many studies on premixed-prevapori zed combustion have been 
made in an attempt to understand better this type of process 
and its potential for applications to future aircraft gas 
turbines. The lean premixing-prevaporizing technique has 
exhibited low emissions in flame-tube studies. For this 
reason, design information and specific studies of various 
aspects of lean, premixed, prevaporized combustion are being 
conducted in a flame-tube experiment to evaluate their 
applicability to aircraft engines. 

The test objective is composed of two parts. The first part 
involves the study of the effects of fuel-air preparation 
characteristics on combustor performance and emissions at . 
temperature and pressure ranges representative of actual gas 
turbine combustors. The second part will investigate the 
effect of flameholding devices on the vaporization process and 
NO^ formation. Flameholder blockage and geometry are some of 
the elements that affect the recirculation zone characteristics 
and subsequently alter combustion stability, emissions and per- 
formance. 

A water cooled combustor is used as the test rig. Pre- 
heated air and Jet A fuel are mixed at the entrance of the 
apparatus. The fuel is injected in two zones. A fully 
vaporized fuel injector is positioned on the centerline of the 
test section and a second non-vapori zed injector is placed 10 
inches from the first one. A vaporization probe is used to 
determine percentage of vaporization and a gas sample probe to 
determine concentration of emissions in the exhaust gases. The 
experimental design will be presented and experimental expected 
results will be discussed. 
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NASA LEWIS FUEL VAPORIZATION STUDY 


OBJECTIVES 


STUDY THE EFFECTS OF FUEL-AIR PREPARATION CHARACTERISTICS ON COOBUSTOR PERFORMANCE 
AND EMISSIONS AT HIGH PRESSURES. 

- DEGREE OF VAPORIZATION 

- DROP SIZE DISTRIBUTION 

- PREMIX QUALITY (FUEL-AIR DISTRIBUTION) 

INVESTIGATE THE EFFECT OF FLAMEKOLDING DEVICES ON THE VAPORIZATION PROCESS AND 
ON NO FORMATION. 

- CHARACTERISTICS OF BLOCKAGE 

- LIQUID COLLECTION 

- REATOMIZATION EFFICIENCY 


EXPECTED RESULTS 


0 THE OBTAINED DATA WILL PROVIDE A BETTER UNDERSTANDING OF THE FOLLOWING SUBJECTS. 

- THE IMPACT OF THE DEGREE OF FUEL VAPORIZATION UPON EMISSIONS AT 
HIGH PRESSURE AND TEMPERATURE. 

- COLLECTION AND REATOMIZATION CHARACTERISTICS OF FLAMEHOLDER DEVICES. 


• MEASUREMENTS HILL BE USED FOR COMPARISON AND DEVELOPMENT OF THEORIES FOR THE 
PURPOSE OF MODELLING THE COMBUSTOR PHENOMENA. 
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EFFECT OF LIQUID DROPLETS ON TURBULENCE 
STRUCTURE IN A ROUND GASEOUS JET 


S.E. Elghobashi and A.M. Mostafa 
Mechanical Engineering Department 
University of California 
Irvine, California 92717 


Abstract 


Experiments [1-5] show that the addition of solid partic.es or liquid droplets 
to free turbulent jets modifies their turbulence structure. For example the 
turbulence intensity and shear stress are reduced due to the presence of the 
dispersed phase. Although the nature of interaction between the dispersed 
phase and the carrier fluid is rather complex and not well understood at present 
[6], there have been several attempts to predict the behavior of particle-laden 
turbulent flows [1,7,8]. 

The objective of the present work is to develop and validate a second-order 
model which predicts the modulation of turbulence in jets laden with uniform 
size solid particles or liquid droplets. 

The approach followed here is to start from the separate momentum and continuity 
equations of each phase and derive two new conservation equations. The first 
is for the carrier fluid's kinetic energy of turbulence and the second for the 
dissipation rate of that energy. Closure of the set of transport equations is 
achieved by modeling the turbulence correlations up to a third order. 

The coefficients (or constants) appearing in the modeled equations are then 
evaluated by comparing the predictions with LDA-measureraents obtained recently 
in a turbulent jet laden with 2-OOu solid particles. This set of constants is 
then used to predict the same jet flow but laden with 50 n solid particles. The 
agreement with the measurement in this case is very good . 
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ASSUMPTIONS 

t BOTH PHASES BEHAVE MACROSCOPICALLY 
AS A CONTINUUM , BUT ONLY THE CARRIER 
FLUID BEHAVES MICRORCOPiCALLY AS A 
CONTINUUM 

• THE DISPERSED PHASE CONSISTS OF 
SPHERICAL PARJICLES OF UNIFORM SIZE' 

• NO COLLISIONS OCCUR BETWEEN 
THE PARTICLES 

t NO PHASE CHANGE TAKES PLACE 
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GCA/ERNING EQUATIONS 
MOMENTUM OF CARRIER FLUID 


^ ^ m - (X - K$2) P 


- F({>2(u^ - v^)] + j ^ j 
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+ PiU^<J>iu^ + * 


MOMENTUM 
OF DISPERSED PHASE 
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+ 8j^*2(P2 - Pi) + f • 


CONTINUITY 


The mean continuity e<juation of the dispersed phase is 


♦z,t + ^ + (<fr2V^) i ■ 0 ; 


the mean global continuity is 
♦ l + <^2 “ 1. 
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TURBULENCE KINETIC ENERGY 
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A SAMPLE CALCULATION 



MEAN VELOCITY 
OF CARRIER FLUID 



0.00 0.05 0.10 0.15 


r/x 

Experiment Predictions 
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Two Phose A 


59 





TWO PHASE FLOW EXPERIMENT 





Vehdt^ Tro^ilfS r/D=:2.0 


6 




0.20 




SPREADING RATE 



Single Phase o 

Two Phase ^ 


CONCLUDING RE.¥.RKS 

• IT HAS BEEN SHOWN THAT THE NEWLY- DEVELOPED TURBULENCE 
MODEL PREDICTS THE BEHAVIOR OF THE TWO- PHASE JET 

• A NUMBER OF UNCERTAINTIES REMAIN AND ONE SHOULD.FOR THE 
PRESENT . EXERCISE DUE CAUTION IN APPLYING THE MODEL TO 
FLOWS THAT ARE SUBSTANTIALLY DIFFERENT FROM THE ONE 
COSIDERED HERE. 

• EXPERIMENT IS URGENTLY NEEDED TO MEASURE MEAN AND TURBULENT 
QUANTITIES IN A GASEOUS JET IN WHICH UNIFORM SIZE DROPLETS 
ARE INJECTED. 
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DILUTION ZONE MIXING 


James D. Holdeman 
NASA Lewis Research Center 
Cleveland, Ohio 44135 


Motivated by considerations of dilution zone mixing in gas 
turbine combustion chambers, NASA sponsored, in 1972 — 1975 
contract and grant studies of the mixing characteristics of a 
row of jets injected normally into a duct flow of a different 
temperature (references 1, 2, & 3) - Based on the favorable 
response to these studies, and the areas for further work which 
they identified, NASA Lewis is currently conducting a balanced 
program of contract, grant, and in-house research on various 
aspects of the 'jet in a confined crossflow' problem. Included 
in these are: 1) development of interactive computer codes for 

analysis of dilution jet mixing, and 2) extension of the 
experiments on jets in a confined crossflow. 


From the data of References 1 Zi 2, an empirical model was 
developed (refs. 4 Sc 5) to describe the observed temperature 
distributions. The current interactive code provides a 3— D 
pictorial representation of the temperature, as given by these 
correlations, for any user— specif ied downstream location, flow, 
and orifice parameters. Although calculations can be performed 
for (almost) any flow and geometric conditions of interest, they 
are, of course, most reliable for conditions within the range of 
the experiments. These codes will be improved and extended, and 
options added, as new data become available. 


The experiments in References 1 to 3 dealt primarily with a 
single row of jets mixing into an isothermal flow in a constant 
cross-section duct. Variations in the mixing were observed as a 
function of jet— to— mainstream momentum ratio, orifice size, and 
spacing. The current experiments examine perturbations of this 
problem characteristic to gas turbine combustion chambers, 
namely: flow area convergence, non— i sothermal mainstream flow, 
and opposed in-line and staggered injection. 


Specifics of grant studies on free— stream turbulence effects and 
reverse flow geometries, and analytical calculations of jets in 
crossflow are presented in companion papers at this conference. 
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B .5 1 

0=tTn-T3/CTn-TJ3 


5/D=Z H/'l>=4 J=32 l>R=2 Cl>=.BH 

X/H=.5 MR=.5 TS=.33 PF=.3H 

Typical Temperature Profile Distribution in Y - 2 Plane 
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x/ll = . 5 ; (S/D)(H/D) = 16 


16 32 



Varations in Temperature Distributions with Momentum Ratio 
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0 1 0 1 0 1 
THETA = (T«, - t)/(T» - Tj) 


Variation in Temperature Distributions with Ori-Fice Size 
at Constant Spacing to Height Ratio (S/H => .5; J — 32> 
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a) Hot Jets — Ambient Mainstream; DR=.75, J-31, R=6.4 



b) Ambient Jets — Hot Mainstream; .DR=:2.2, J=26, R=3.5 



c> Ambient Jets - Hot Mainstream; DR=2.3, J=109, R=6.9 
Effect of Density Ratio on Temperature Profiles (S/D=4, H/D=B) 
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a) Hot Jets - Ambient Mainstream; DR=.62, -J=31, M^S.S 



b) Ambient Jets — Hot Mainstream; DRs2.1, J=:21.6, M=»6.3 



c) Ambient Jets - Hot Mainstream; DR=2.1, J=5.7, M=3,2 
E-f^ect of Density Ratio on Temperature Profiles (S/D*=2, H/D=4) 
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Infl 


X/HB « 0.600 
X/DJ 3 2.68 


X/-H0 3 0.750 
X/OJ 3 3.87 


X/m 3 1.000 
X/OJ 3 6.18 




a) Straight- Duct 



b) Symmetric Convergence (.5/1) 



c) Injection Wall Converging (.5/1) 

Influence of Flow Area Convergence on Temperature Profiles 


5 



a) Single— side (*top) Injection <S/H— *5, J— 26.3) 



iimohfi/tratii-TJi in«i»-ti/tnfiDi>rjt 


b) Two-side (opposed) Injection (S/H=.25, J=25) 
Comparison Between Single— side and Opposed Jet Injection 
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INFLUENCE OF A LARGE FREE STREAM DISTURBANCE LEVEL ON 
DYNAMICS OF A JET IN A CROSS FLOW 


John J. Foss and Candace E. Wark 
Department of Mechanical Engineering 
Michigan State University 
East Lansing, Michigan 48824 


INTRODUCTION 

Jets of relatively high speed (V.) and low temperature (T.) are 

J J 

often used, in a gas turbine combustor, to cool the gases and to quench 

the chemical reactions. The prior studies, of this particular class 

of the "jet in a cross flow" problem have commonly been executed 

2 2 

with the specific momentum flux ratio [PjVj /PqU^ ] as the primary 
variable, values of 10-100 characterize the range of technological 
interest for the combustor cooling problem. 

The jet in a cross flow may be subdivided into two general 
regions: "interacting" and "downstream". Our interest is in the 
former and in the physical agents that are responsible for the: "jet 
turning into the streamwise direction" and the mixing of the jet and 
the cross stream fluid. A representation of the interaction region 
for isothermal mixing at (V-/U ) = 3.0 is shown in Figure 1; this 

vJ 

figure, from Foss [1980], serves to characterize the interaction 
region. 

The prior research studies of this problem have used wind tunnels 
for the cross stream; a concomitant attribute of these studies has been 
the presence of a low free stream disturbance level. The purpose of our 
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investigation is to provide direct observations of the jet trajectory 
and mixing in the presence of a large disturbance condition. There 
are two distinctive features of our investigation. These are described 
below. 

The flow facility, of the MSU Free Shear Flows Laboratory 
(Figure 2), provides a large, planar, shear layer for use in the present 
study. Specifically, we have placed the jet (d. = 10mm) such that it 
exhausts into the middle region (u/U^ = 0.5) of the shear layer at the 
end of the 3 meter test section. The local vorticity thickness (6 ) is 
large with respect to the jet diameter - 58. Based on this jet 

diameter, the turbulence filed is essentially homogeneous. (The 
gradients of the mean velocity and the turbulence intensity are con- 
sidered to cause second order effects; the influence of these dis- 
tributions could be examined with appropriate numerical models of the 
flow.) The details of the experimental facility are shown in Figures 
3 to 4 . 

The second distinctive feature of our study is the use of an 
array of 76, fast response (t= 0.1 msec) thermocouples to document 
the instantaneous temperature field at the "end" of the interaction 
region. Specifically, the array will be placed in the flow, as shown 
in Figure 3 and the complete set of simultaneously sampled thermo- 
couple voltages will provide discrete values to characterize the 
temperature field: T(y, z, t.).. The temperature field measurements 

will be repeated such 2.5 msec; hence, both the instantaneous field 
and its temporal evolution will be evaluated. Fpr comparison, the 
experiments will be repeated in the undisturbed, high speed, region of 
the flow field. 
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The "planes" of the instantaneous temperature field provide a 
unique data base for the evaluation of the trajectory and the mixing 
of a "jet in a cross flow." The temperature data, in the undisturbed 
cross stream, will be used to extend our understanding of this basic 
flow field. The comparative measures: with and without a disturbed 
cross stream, will be used to identify the influence of this dis- 
turbance parameter on this flow field. 
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Figure 1. Schematic representation of a **large R** jet in a cross 
flow from Foss [1979] • 

Notes: schematic representation of the observed shear layer 

instability 

schematic representation of the vorticity vector. 

Surface stress and A,B,C,d fluid trajectory lines have been traced 
from the appropriate photographs. . .their positions are shown to 
proper scale with respect to the jet hole diameter. 

The “shroud" of fluid, which covers the flanks of the wake region 
and which is marked to the sheared fluid from the A and B inputs, 
is not shown. 

The forward stagnation nodes (N 2 ,N 2 ), as shown in Figure 33, exist 
for the large R condition, they are only shown explicitly in that 
figure. 

The 3l> CX2» 32 streamlines and the cross hatched section of the 
jet represent conjectural estimates of the flow behavior. The forme 
are to show the presence of a sharp division between jet and cross 
stream fluid and the formation of the bound vortex from the inter- 
action of these two streams. 
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Figure 3. Schematic representation of the flow field, the 
Thermocouple Sampler, and Instrumentation 
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DILUTION JET EXPERIMENTS IN 
COMPACT COMBUSTOR CONFIGURATIONS 


Isaac Greber and James Zizelman 
Department of Mechanical and Aerospace Engineering 
Case Western Reserve University 
Cleveland, Ohio 44106 


This project concerns the effects of cooling jets on the velocity and tem- 
perature fields in a compact reverse flow combustor. The work is motivated by 
the need to limit the temperatures of past-combustion gases in jet engines to 
values within the endurance capabilities of turbine blades. The application 
requires not only that the temperature be kept sufficiently low but also that 
a suitably tailored temperature profile be provided at the combustor exit, with 
higher temperatures generally permissible at the blade tip than at the blade 
root because of higher centrifugal loads at the root. 

Flows in reverse flow combustors accelerate both longitudinally because of 
area changes and transversely because of flow turning. The current project 
started with flow visualization experiments in water, using aqueous solutions 
of zinc bromide to model the relatively higher density of cooling jets. These 
flow visualization experiments were conducted in simple two dimensional con- 
figurations designed to examine separately the effects of longitudinal and 
transverse acceleration. The next phase of the work consisted of temperature 
measurements in a model reverse flow combustor, using a rake of thermocouples 
which could be moved both longitudinally and transversely to provide transverse 
temperature profiles at a number of longitudinal stations. Single jets, and 
rows of cooling jets of different spacing, were injected at several locations 
with varying flow rates and temperatures, to produce a useful range of momentum 
ratios and density ratios. In addition, a semi-empirical calculational model 
was developed to predict the behavior of a single cooling jet in the reverse 
flow combustor configuration. 

Results to date show that single jet temperature trajectories are swept 
toward the inner wall of the turn, whether injection is from the inner or outer 
wall. A widely spaced row of jets produces a trajectory similar to that of 
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a single jet. As spacing is reduced, jet penetration is also reduced, and the 
cooling jets tend to remain close to the wall from which they are injected. 

The results to date suggest that suitable cooling and temperature distribution 
tailoring can be accomplished without injecting cooling jets upstream of the 
turn, and thus it appears that combustors can be made significantly smaller 
than current designs. 

The current phase of this work is directed at the acquisition of velocity 
measurements as well as temperature measurements in the combustor. To this 
end, a combined pitot-static tube and thermocouple probe has been constructed 
and tested, and a rake of such tubes has been designed and is being constructed. 

In order to keep new temperature measurements to a minimum, testing will be 
carried out in such a way as to attempt to reproduce earlier temperature conditions. 

An attempt to develop a semi-empirical calculational model for a row of 
jets, similar to that for a single jet, was limited by the absence of fundamental 
mixing information for rows of jets. Consequently, experiments were performed 
to obtain entrainment rates for rows of jets; the results are currently being 
analyzed. Using these results, an attempt will be made to complete the calcu- 
lational model. 

Throughout this work we have been fortunate to have the support and advice 
of Steve Riddlebaugh of the NASA Lewis Research Center, and we are pleased here 
to thank him. 


PUBLICATIONS 

Lipshitz, A. and Greber, I., "Turbulent Jet Patterns in Accelerating Flows," 
AIAA Paper 81-0348, 1982. 

Lipshitz, A., "Dilution Jets in Accelerated Cross Flows," CWRU, 1981. 

Riddlebaugh, S., Lipshitz, A. and Greber, I., "Dilution Jet Behavior in the 
Turn Section of a Reverse Flow Combustor," AIAA Paper 82-0192, 1982. 
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Figure 3. - Cutaway sketch of test combustor. 



measuring stations. 
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PROFILES 

Figure 5 - Temperature distribution in the combustor 
turn without dilution jets. 



DENSITY RATIO 2.74 
MOMENTUM RATIO 
A 5.69 

a 7.67 
o 11.7 


Figure 6. - Jet trajectories showing effect of momentum 
ratio on a jet injected from the inner wall. 
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DENSITY RATIO 2.75 
MOMENTUM RATIO 
5.67 

□ 7.65 
o 9.81 


Figure 7. - Jet trajectories showing effect of momentum 
ratio on jet injected from the outer wall. 
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— o— MEASURED 
PREDICTED 


Figure 8. - Comparison with model for a single jet injected 
from the inner wall. 


90 



10 



MEASURED 

PREDICTED 


Figure 9. - Comparison with model for a single jet injected 
from the outer wall . 
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Figure 10. - Behavior of a widely spaced row of jets 
(spacing ratio of 7.41) injected from the inside wall. 


91 






^ ,< _ . ^ 




EXHAUST 


STATION 5 


SINGLE JET 


JET 

INJECTION 

POINT 


STATION 3 


COMBUSTION 
GASES m 


STATION 1 

RADIAL TEMPERATURE 
PROFILES 


JET CENTERLINE TRAJECTORY 


Figure 11. - Behavior of a closely spaced row of jets 
(spacing ratio 2.47) injected from the inside wall. 
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ANALYTICAL CALCULATION OF A SINGLE JET 
IN 

CROSS- FLOW AND COMPARISON WITH EXPERIMENT 


R. W. Claus 

NASA Lewis Research Center 
Cleveland, Ohio 44135 


INTRODUCTION 

With the increasing costs of combustor development testing, 
a great deal of interest has focused on the use of numerical 
models to screen deign changes or develop new combustor 
concepts. This type of combustor based design process now 
appears to be feasible through the use .of 3D combustor 
performance models such as those reported in refs. 1 and 2. 
Ultimately, actual hardware testing may, someday, be used only 
for final design verification. 

A number of major restrictions must be overcome before this 
type of design methodology can be adopted. First, the proper 
physics must be incorporated into the differential equations 
used in the combustor model. Second, numerical methods must 
accurately solve these differential equations. Finally, the 
accuracy of the resulting code must be assessed against 
fundamental data and improvements made to the code to alleviate 
identified deficiencies. 

Currently available 3D combustor performance models have 
yet to be thoroughly assessed. A few comparisons have been 
made against actual combustor hardware, ref. 3, but these have 
not been conclusive. A more logical first step is to examine 
the extent to which three-dimensional hydrodynamic processes 
can be calculated. One flow field of this type for which a 
great deal of experimental data exists is jets in cross-flow. 

Jets in cross-flow are practically revelant to the gas 
turbine combustor designer. Cooling air jets (dilution jets) 
are used to control the hot gas temperature profile entering 
the turbine. They are also used to set up aerodynamic patterns 
within the combustor which promote mixing and control local 
burning zone stoichiometry. As a result the jet penetration 
and mixing characteristics of jets in cross-flow are of primary 
concern in the combustor design process. 

There have been a number of previous calculations of jets 
in cross-flow, refs. 4, 5, and 6, however, these studies have 
been limited by two main factors. First, although a great deal 
of experimental data exists, rarely have important parameters 
such as the turbulence field, inlet velocity profiles and jet 
mixing characteristics been fully measured. This limits the 
flow field quantities one can compare and imposes the need to 
assume inlet boundary conditions of the calculation. The 
recent measurements of reference 7 greatly reduce this pro- 
blem. Second, core storage and economy requirements have 
limited previous calculations to coarse grid systems. This 
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results in some numerical error being present in the computed 
solution, which can possibly call into question any conclusions 
drawn from these, studies. 

The present report expands on this previous work by 
employing a series of progessively finer grid systems to 
calculate the single jet in cross-flow experimentally measured 
in ref. 7. These experimental measurements provide a fairly 
complete collection of velocities, turbulence intensities, and 
jet concentration profiles with measurements of the inlet 
field. The use of a series of progressively finer grid systems 
allows a differentiation between numerical errors and the 
hydrodynamic modeling assumptions embodied in the 3D combustor 
code of ref. 1. The results of this comparison will provide 
additional insight into the def iiciencies of the turbulence 
model and the code numerics. 


CONCLUSIONS 

Employing a 3D finite-difference model to analyze the jet 
flow field of ref. 7, the following conclusions were determined. 

1. With a reasonable number of grid points (approximately 
40 X 30 X 20) calculated jet penetration and con- 
centration profiles agreed well with experimental 
measurements . 

2. For the cross-stream vortex, the agreement between 
experimental and calculated results become less 
qualitative as additional grid points were added - 
indicating a deficiency of the isotropic turbulence 
model . 

3. The calculated results of the finest grid examined 
(90 X 40 X 22) were grid dependent. An improved 
numerical scheme is required to remove the effects of 
numerical diffusion for the flow geometry examined. 
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NUMERICAL MODELING 




COMBUSTION 

FUNDAMENTALS 


OBJECTIVE: ASSESS AND IMPROVE THE STATE OF THE ART 

IN COMBUSTION MODELING 


APPROACH: EMPLOY FINITE DIFFERENCE MODELS OF THE 

TIME-AVERAGED NAVIER STOKES EQUATIONS 


BENEFIT: AN INCREASED DESIGN CAPABILITY WITH 

REDUCED DEVELOPMENT COSTS 
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DILUTION JETS - 3D FLOWFIELD 


COMPUTATIONAL DETAILS: 

20x20x12 grid 20 CPU minutes 

40x30x20 grid 2 CPU hours 

90x40x22 grid 10 CPU hours 

RESULTS : 

QUALITATIVE PREDICTION OF 

1. JET PENETRATION 

2. MIXING CHARACTERISTICS 

LIMITATIONS INCLUDE 

1. GRID DEPENDENT SOLUTIONS 

2. ISOTROPIC TURBULENCE MODEL 


y 


A 



z 


THREE DIMENSIONAL SINGLE, FREE JET FLOW FIELD SCHEMATIC 
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AXIAL VELOCITY PROFILES ALONG JET CENTERLINE AT X/D=2 



CALCULATED AND EXPERIMENTAL lET TLUIO CONCENTRATION PROFILES AT X/D<8 


20x20x12 



40x30x20 

EXPERIMENT CALCULATION 



90x40x22 

EXPERIMENT CALCULATION 
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Y-Z PLANE VELOCITY VECTOR PLOTS OF TWO DIFFERENT CALCULATIONS 
OF THE SINGLE JET FLOW FIELD 


20x20x12 


AT APPROXIMATELY X/D- 0.7 

40x30x20 


VJD 


CALCULATION 

r 

vnu\/ULniiuM 




• • * \ 1 



- " M 



\ 



\ \ t 


. . . : ^ ' 
* » • ' * 


^ ill • 4 -^'1 

M : 


0 3 

ZiD 



CS-62-34S3 


40x30x20 GRID CALCULATION 

AXIAL VELOCITY PROFILES ALONG JET CENTERLINE TURBULENCE AND JET FLUID CONCENTRATION 
ATX/d ■ 6 PROFILES ALONG JET CENTERUNE ATX/D • 6 
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COMPARISON OF TERMS IN THE AXIAL MOMENTUM EQUATION 


X/D-L3 

CONVECTION 


PRESSURE GRADIENT (PG) 

CONVECTION TRl^CATION 




TERMfPG^, Uflt. 

CS-82-3433 
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TURBULENT COMBUSTOR FLOWFIELD INVESTIGATION 


David G. Lilley 

School of Mechanical and Aerospace Engineering 
Oklahoma State University 
Stillwater, Oklahoma 74078 


Experimental and theoretical research is being undertaken on 2-0 
axisymmetric geometries under low speed, nonreacting, turbulent, swirling flow 
conditions. The flow enters the test section and proceeds into a larger 
chamber (the expansion ratio D/d = 2) via a sudden or gradual expansion (side- 
wall angle a = 90 and 45 degrees). Inlet swirl vanes are adjustable to a 
variety of vane angles with values of = 0, 38, 45, 60 and 70 degrees being 
emphasized. The objective is to determine the effect of these parameters on 
isothermal flowfield patterns, time-mean velocities and turbulence quantities, 
and to establish an improved simulation in the form of a computer prediction 
code equipped with a suitable turbulence model. This is a prerequisite to the 
prediction of more complex turbulent reacting flows, and successful outcomes 
of the work can be incorporated into more combustion- and hardware-oriented 
activities, including incorporating the modeling aspects into already existing 
comprehensive numerical solution procedures. 

New features of the present year's study are; the turbulence 
measurements are being performed on swirling as well as nonswirling flow; and 
all measurements and computations are also being performed on a confined jet 
flowfield with realistic downstream blockage. Recent activity in the research 
program falls into three categories: 

1. Time-mean flowfield characterization by five-hole pitot probe 
measurements and by flow visualization. 

2. Turbulence measurements by a variety of single- and multi-wire hot- 
wire probe techniques. 

3. Flowfield computations using the computer code developed during the 
previous year's research program. 

In the experimental approach, flow visualization is accomplished by still 
and movie photography of helium-filled soap bubbles and smoke produced by an 
injector and a smoke wire.^ Time-mean velocities have been measured with a 
five-hole pitot probe at low swirl strengths, ^ and this has recently been 
extended to higher swirl strengths, including the effect of size and location 
of a downstrean blockage of area ratio 2 and 4. ^ Turbulence measurements have 
recently been completed on swirling (up to (J) = 45 degrees) as well 
nonswirling flows using a six-orientation single-wire hot-wire technique. 
Although all Reynolds stresses are deducible by this technique, a new cross- 
wire probe is currently being used for more accurate shear stress 
evaluation. A series of hot-wire measurements are also in progress, using a 
three-wire hot-wire probe with direct computer inteface and data reduction. 
These data will enable turbulent viscosity components to be deduced (from the 
turbulence stress vs. mean velocity spatial gradients) and will aid the 
turbulence modeling aspects of the study. 
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In the computational approach, an advanced computer code has been 
developed to predict corresponding confined swirling flows to these studied 
experimentally.® Tentative predictions were given earlier.^ Predictions have 
now been made which include realistic inlet conditions for a complete range of 
swirl strengths. Predictions of downstream blockage (two blockage sizes at 
two axial locations) have also been completed and these are being assessed in 
the light of the experimental data. ® Advanced versions of the code will 
incorporate modification in the form of state-of-the-art turbulence models, 
deduced from the experiments and elsewhere, so as to accommodate better the 
physics of swirling recirculating flows. 

Proposed future studies include: expansion ratio D/d = 1.5, 
correspondingly predictions using advanced turbulence models, inlet turbulence 
level effects to be clarified, turbulence length scales to be measured, and 
turbulence model developments to be made. The completion of these tasks will 
make a substantial contribution to the understanding and predictive capability 
of complex turbulent swirling recirculating flows. 
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1 . 2 Research Objectives 

1 . To dete rmi ne the effect of 

* swi rl vane angle (}> 

* side-wall angle a 

* downstream blockage area ratio AR 

and location L/D 
on : 

* isothermal flowfield patterns 

* time-mean velocities 

* turbulence quantities 

2. To establish an improved simulation in 
the form of: 

* computer prediction code 

* suitable turbulence model 


1 . 3 Research Approach 


1. Time-mean flowfield characterization by five- 
hole pitot probe measurements and by flow 

vi sual i zati on . 

2. Turbulence measurements by a variety of single- 
and multi-wire hot-wire probe techniques. 

3. Flowfield computations using the computer code 
developed during the previous year's research 
program. 


1.4 Progress During First Year (1980-1981) 

1. Test f aci 1 i ty - design and construction, including 

variable angle swirler and expansion blocks. 

2. Experimental techniques - flow visualization, five- 

hole pitot probe and reduction computer code. 

3. Computational code development - advanced version 

of a pri mi t i ve- va ri ab 1 e solution procedure - 
STARPIC. 

4. Flowfield characterization - emphasis on time-mean 

flow character, using flow visualization, five- 
hole pitot probe and tentative predictions. 
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1.5 Progress During Second Year (1981 - 1982) 

1. Higher swi r1 s trengths - flow visualization and 
five-hole pitot probe measurements. 

2. Swirler effectiveness - exit profiles of u, v, w, 
and p accurately established. 

3. Flowfield predictions - the need for specifying the 
inlet conditions very precisely. 

4. Downstream blockage - effects of area ratio and 
axial location studied by flow visualization, five- 
hole pitot probe and corresponding computer 

predi cti ons . 

5. Turbulence measurements - in nonswirling and swirling 
flow via one-wire, two-wire and three-wire hot-wire 
methods . 


1.6 Present Studies (1982 - 1983) 

1. Expansion ratio D/d = 1.5 

2. Corresponding predictions 

3. Turbulence measurements at higher swirl 
strengths, with and without downstream blockage. 

4. Turbulence length scale being measured 

5. Turbulence model developments 
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3. Time-Mean Velocity Measurements 
3.1 Swirl Effect on Axial and Swirl 


x/D 


0 1 





(c) (!) = 70° 




Ve 1 oci ty Fields 





f 
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3 . 5 Swirl Effect on Streamlines 
Blockage at L/D = 2.0 


Weak Downstream 





Swirl Effect on Strea 
Blockage at L/D = 2.0 


ines; Strong Contraction 













6 . Cl os u re 

1. Measurements of effect of swirl, expansion angle 
and contraction blockage ££ flowfield patterns, 
time-mean velocities and turbulence quantities. 

2. Computer prediction studies of associated phenomena. 

3. Turbulence model developments for the simulation 
of swirling recirculating flow. 


117 




SMALL GAS TURBINE COMBUSTOR PRIMARY ZONE STUDY 


R. E. Sullivan and R. D, Sutton 
Detroit Diesel Allison 
Ind ianapol is , Indiana 4 62 02 

SUMMARY OF PRESENTATION 

This talk will summarize the combustion research program sponsored by 
NASA Lewis Contract No. NAS3-22762, "Small Gas Turbine Combustor Primary 
Zone Study." Recent publication on this effort was the June 22, 1982 
AIAA Paper 82-1159 authored by R. Sullivan, A. Novick, G. Miles of 
Detroit Diesel Allison, and by D. Briehl of NASA Lewis Research Center, 
who is the Project Manager. 

The presentation will describe the basic elements of a design methodology 
program to obtain the maximum performance potential of small reverse- 
flow annular combustors. Three preferred combustion design approaches 
for internal flame stabilization patterns were selected for study. 

Design features were incorporated in the combustors to address the 
performance limiting problem areas associated with small annular 
combustors. Performance was predicted using a 3-D aerodynamic/chemical 
kinetic elliptic flow analysis, initially developed by Garrett Corporation 
for the USARTL- The analytical performance predictions are compared 
with actual test results, measured at the exit plane of the primary 
zone. The findings illustrate that the analytical flowfield predictive 
models provide a very useful design tool for understanding the combustion 
performance of a small reverse-flow annular combustor. 
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OBJECTIVE 


• FORMULATE UMDERSTANDING OF PRIMARY ZONE 
AERODYNAMICS 

• RELATE TO PERFORMANCE OPTIMIZATION 

• mPROVri DESIGN METHODOLOGY OF REVERSE FLOW 
ANNULAR COMBUSTORS 

• INTERACTIVE DESIGN. ANALYSIS AND TEST 


CHARACTERISTICS OF 
REVERSE FLOW ANNULAR COMBUSTOR 



CONCERNS: 

• COOLING 

• ATOMIZATION 

• FUEL IMPINGEMENT 

• WALL OUENCHING 

• REACTION TIME 
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COMBUSTOR PRIMARY ZONE CONCEPTS 


DOUBLE VORTEX REVERSE COUNTERFLOW 



SINGLE VORTEX 



COMBUSTOR EXPERIMENTAL MEASUREMENTS 



ALL CONCEPTS 


COMBUSTOR OUTLET 
360 TRAVERSE PROBES 
PRESSURE PROBE 
4 ELEMENT -I 




PRIMARY ZONE 
PROBE LOCATIONS 
4 PROBES WITH 
4 RADIAL PORTS 


TEMPERATURE RAKE 
4 ELEMENT 

2 PLACES, 180“APART 


GAS SAMPLING PROBE 
3 RADIAL DEPTHS 


3 
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EMISSION INDEX-CH. COMBUSTION EFRCIENCY-% 


TYPICAL TEST RESULTS 
CONCEPT I COMBUSTOR 
EFFKIBICY AND EMSSHNS 






PRIMARY ZONE PROBE LOCATIONS 
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EMISSION INDEX-NO , EMISSION INDEX-CO 


TYPICAL TEST RESULTS 
CONCEPT I COMBUSTOR 

PATTERN FACTOR 




-I 1 I I 

25 50 75 100 

POWER-V. TI82-3421 


TYPICAL TEST EVALUATION AND 
CORRELATION WITH ANALYSIS 
CONCEPT I COMBUSTOR 
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80 
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EXHAUST (REF) 

^>-^PZ PROBES 


NOZZLE MID BETWEEN 

CENTERLINE LOCATION NOZZLES 

PROBE LOCATIONS 


: 12 
10 
8h 
6 
4 
2 
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EXHAUST (REF) 

>PZ PROBES 


NOZZLE MID BETWEEN 

CENTERLINE LOCATION NOZZLES 

PROBE LOCATIONS 



u 20 h 


CO g [ - i » 

NOZZLE MID BETWEEN 

CENTERLINE LOCATION NOZZLES 

PROBE LOCATIONS 


PREDICTED MARC-I 



CENTERLINE LOCATION NOZZLES 
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ANALYTICAL DESIGN PROCEDURE 


• INCORPORATE ANALYTICAL ANALYSIS TO PREDICT 
RIEL-AIR DISTRIBUTION 

• DERNE PRIMARY ZONE AERODYNAMICS 

• TAILOR RIEL PLACEMENT TO AIR PATTERNS 

• CORRELATE TEST RESULTS TO UPDATE DESIGN 
PROCEDURES 

• USE ANALYTICAL ANALYSIS TO GUIDE 
COMBUSTOR M0DIRCAT10NS 

• RELATE OVERALL PERRIRMANCE TO 
CONDITIONS IN PRIMARY ZONE 


MARC-I CODE DESCRIPTION 

• 3-D AERODYNAMIC/REACTING ELLIPTIC FLOW ANALYSIS 

• ARMY/GARREH CODE USED AS BASELINE 

• PRIMITIVE-VARIABLE. FINITE-DIFFERENCE CODE 

• SOLVES NAVIER-STOKES EQUATIONS IN 3-D 

MARC-I CODE ELEMENTS 

• K-6TURBULENCE MODEL 

• FUEL SPRAY VAPORIZATION MODEL 

• TWO-STEP REACTION MODEL BASED ON ARRHENIUS AND 
EDDY-BREAKUP CONCEPTS 

Cx^Y * CO '*C02 
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TYPICAL COMPUTER OPERATION 



« ^ ^ ^ 
y» / / 




>• 

>• 









lONemiDINU CROSS-SECTION 


SECTOR CROSS-SECnON 


GRID SIZE; 17X13X13 

REQUIRED ITERATIDHS: 3D0 (NEW DESIGN) 

IDO (DESIGN REVISION) 

COMPUTER TIME (IBM 370): 15 MINUTES ($100/100 ITERATIONS) 
GRID SPACING: aiS- 020 INCH 


ODA M0DIRCAT10NS 

• GEOMETRIC REHNEMENT 

• PRECHAMBERS 

• SWIRLERS 

• VARIABLE OOME SHAPES 

• COOUNG AIR ADMISSION 

• IMPROVED DATA PRESENTATION 

• VELOCITY VECTOR VISUAUZATION 

• 3-0 CONTOUR puns 

• RESTART CAPABILITY 

• DATA STORAGE FOR SIMILAR DESIGNS 
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ANALYTICAL TECHNIQUE 

• MARC-I 

• VELOCITY VECTOR AND TEMPERATURE CONTOUR PLOTS 



(B) TEMPEBATURE CONTOUR PLOT 


MARC-I 

• TEMPERATURE CONTOURS (DEG F/lOO) 



VELOCITY VECTOR PLOTS PNCEFT II) 




-=t 



|A| LONGUDDiNAL SECTION THROUGH AXIAL SWIRL£R AND PRIMARY HOLES 



(B| LONGITUDINAL SECTION THROUGH INTERMEDIATE HOLES 


VELOCITY VECTOR PLOTS (CONCEPT IIQ 




127 



FUEL-AIR RATIO 


PERFORMANCE OPTIMIZATION OF 
COMBUSTOR MODIHCATIONS 

• ANALYTKAL ANALYSIS. MARC-L USED TO GUIDE 
COMBUSTOR CONCEPT I MODIFICATIONS 


CONCEPT I • MODIFIED 

PREDICTED 

MEASURED ^ 



CIRCUMFERENTIAL LDCATION- DEGREES 
RADIAL AVERAGING 


CONCEPT I - MODIRED 

PREDICTED 

MEASURED 


22.5 DEG 15.5 DEG 11.25 DEG 

0.16 
^ 0.15 

CO 
2 

S 0.14 

oe 

. 0.13 

0 0.05 0.10 0 0.05 0.10 0 0.05 0.10 

FUEL-AIR RATIO 

RADIAL VARIATION 



CONCEPT I - MOD I - 80-PERCENT POWER 
FUEL-AIR RATIO 



*SAMPUNG REFERENCE 


128 



APPUCATION OF MARC-I TO DESIGN 



IS 



ID) CONCEPT III BASEUNE SECTOR NITN FUEI-AIR PROmE AT PRIMARY ZONE EXIT 




|B) CONCEPT III MODIRCATION (ADDENDUM) 


SUMMARY 

• USED 3 D MODEL TO DEFINE PRIMARY ZONE AERODYNAMIC 
AND FUEL DISTRIBUTION PATTERNS 

• EXPLORED POTENTIAL OF COMBUSTOR CONCEPTS DESIGNED 
TO ADDRESS PROBLEMS OF REVERSE-FLOW ANNULAR 
COMBUSTORS 

• OBTAINED PRIMARY ZONE GAS SAMPLES USING EFFECTIVELY 
DESIGNED WATER-COOLED PROBE 

• CORRELATED PRIMARY ZONE MEASUREMENTS TO ANALYTICAL 
PREDICTIONS AND OVERALL PERFORMANCE 
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MASS AND MOMENTUM TURBULENT TRANSPORT EXPERIMENTS 


B. V, Johnson 

United Technologies Research Center 
East Hartford, Connecticut 06040 


An experimental study of mixing downstream of coaxial jets discharging in an 
expanded duct was conducted (1981) to obtain data for the evaluation and improvement 
of turbulent transport models currently used in a variety of computational procedures 
throughout the propulsion community for combustor flow modeling. The study used 
laser velocimeter (LV) and laser induced fluorescence (LIF) techniques to measure 
velocities and concentration and flow visualization techniques to qualitatively 
determine the time dependent characteristics of the flow and the scale of the 
turbulent structure , 

Qualitative and quantitative studies were conducted of the flow downstream of 
coaxial jets discharging into an expanded duct. The ratio of annular jet diameter 
and duct diameter to the inner jet diameter were .2 and 4, respectively. The inner 
jet peak velocity was approximately one^half the annular jet peak velocity. RCsSuTts 
from the studies were related to the four shear regions within the duct: (1) w„..e 

region downstream of the inlet, (2) shear layer between the jets, (3) recirculation 
region, and (4) reattachment region. 

Flow visualization studies were conducted using dye as a trace material and 
high-speed motion pictures to record the dye patterns in selected r-z and r-f 
planes. The principal results from this study showed: 

1. The larger scales of the turbulent structure were observed to grow from the 
width of the wake region downstream of the inner jet tube to a large fraction of the 
duct diameter immediately upstream of the reattachment zone. 

2. The turbulent eddies were not axisymmetric or periodic at any location 
within the duct. The large scale waves and eddies appeared to have a range of 
wave lengths . 


A detailed map of the velocity, concentration, mass transport rate and momentum 
transport rate distribution within the duct was obtained to provide data for the 
evaluation and improvement of turbulent transport models. Data sets of two velocity 
components pairs were obtained simultaneously to determine momentum transport rate 
and velocities. Data sets of velocity and concentration pairs were obtained 
simultaneously to determine mass transport rate, concentration, and velocity. 
Probability density functions (p.d.f.s) of all the forementioned parameters were 
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obtained from the data sets. Mean quantities, second central moments, skewness and 
kurtosis were calculated to characterize each data set. The principal results from 
this study are: 

3. The axial and radial velocity profiles documented the changes in the shear 
regions within the duct. 

4. The mean and fluctuating concentration profiles documented the inner jet 
fluid distribution within the duct. 

5. The turbulent momentum transport rate measurements in the r-z plane 
documented the local momentum fluxes due to turbulent mixing. Correlation 
coefficients were determined for each measurement location and data set, 

6. Countergradient turbulent axial mass transport was measured in the shear 
region between jets. The peak axial mass transport rates were greater than the peak 
radial mass transport rates even though the axial concentration gradients were 
approximately one-fifth the radial gradients. 

7. The countergradient turbulent axial mass transport was related to the 
general direction of the eddies between the inner and annular jets. The counter- 
gradient axial mass transport occurred when the annular jet was accelerating the 
inner jet fluid. 

8. The skewness and kurtosis of the momentum transport p.d.f.s in the peak 
shear region were approximately the same as previously measured in turbulent 
boundary layers. However, the kurtosis in the low shear region was greater than 
previously measured in the wake region of the turbulent boundary layers. 

9. The peak values of kurtosis for the mass transport p.d.f.s was greater 
than the peak values for the momentum transport p.d.f.s, 

10. The kurtosis for all the transport rate p.d.f.s were an order of magnitude 
greater in the low transport rate regions, including the recirculation region, than 
in the high transport rate regions. 

Report: Mass and Momentum Turbulent Transport Experiments with Confined Coaxial 

Jets by B. V. Johnson, J, C. Bennett, NASA Contractor Report NASA CR-165574, 

November 1981. 

Currently swirling coaxial jets are being experimentally studied using the same 
LV/LIF measurement technique used in the nonswirling coaxial jet study. The only 
change in the flow condition from the previous study is the addition of 30 deg 
mean angle swirl to the annular jet fluid. The flow visualization technique has 
been improved to show interaction of shear layers on both sides of the annular 
stream. 


132 




Figure i. Sketch of inlet and test section; L = 1016 mm, 
R ^1 = 12.5 mm, R-j 2 = 15.3 mm, Ra = 29.5 mm, 

Rq = 61 .0 mm. 


WATER 

BOX 



Figure 2. Optical arrangement for flow 
visualization in r-z plane. 
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RECIRCULATION REGION 



Figure 3. Shear regions of nonswirl inq 
coaxial jets in an enlarged 
duct. 



Figure 4. 


Shear regions of coaxial jets 
with swirl in an enlarged duct. 
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Figure 5. Sketch of optics arrangement. 
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Figure 6. Probability density functions for axial velocity 
and concentration at r/Rg = 0.2; Au = 0.1 m/s, 
u = 1.20 m/s, u' = 0.22 m/s, = -0.15, Ky = 216; 

Af = 0.02, f = 0.28, f = 0.18, Sf = 0.91, Kf = 4.0. 
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Figure 7. Probability density functions for r-z momentum and 
radial mass transport ratio at r/Rp = 0.2; Auv = 
0.01 m2/s2, uv = 0.0142 m^/s2, ctuv = 0.0^ m2/s^, 
Suv = -2.14, Kyy = 11.9; Oyf = 0.04 m/s, uf = 

0.0165 m/s, = 0.028 m/s, Syf = 2.23, Kyf = 10.4 
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Figure 8. Moments of axial velocity profile; o = 
0 deg -OA, iso -•A, 90 270 
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RADIUS RATIO — t/Rq 

Figure 9. Moments of concentration profiles; o = 0 deg -Afck, 
180 -A k 90 - 270 - 0. 



RADIUS RATIO — r/Rg 


Figure 10. Moments of turbulent momentum transport profiles 
0 = 0 deg - O, 180 • 
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Figure 11. Moments of turbulent radial mass transport pro- 
files; 0 = 0 deg - 0 , 180 


138 







m/s 





MASS TRANSPORT 



0 1 2 3 4 5 

z/R 


Figure 13. Regions of countergradient axial m^s transport 

and axial velocity accelerations; uf < 0 is counter- 
gradient region. 
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VELOCITY VISUALIZATION IN GASEOUS FLOWS 


R. K. Hanson, J. C. McDaniel and B. Hiller 
Department of Mechanical Engineering 
Stanford University 
Stanford, California 94305 


Techniques yielding simultaneous, multiple-point measurements of 
velocity in reacting or nonreacting flowfields have the potential to 
significantly impact basic and applied studies of fluid mechanics. This 
research program is aimed at investigating several candidate schemes which 
could provide such measurement capability. The concepts under study 
utilize laser sources (to illuminate a grid, a plane or a volume in the 
flow) which lead to scattered light (from Mie scattering, fluorescence on 
chemiluminescence) which can be monitored by a film-based camera or a 
multi-element solid-state camera (100 x 100 array of photodiodes). Anti- 
cipated experiments will be discussed, and recent results for velocity 
measurements in supersonic flows using a novel Doppler-modulated fluorescence 
concept will be presented. Related research at Stanford on species visuali- 
zation sponsored by AFOSR will also be summarized as time permits. 


OUTLINE OF PRESENTATION 

• MOTIVATION (NEED FOR MULTIPLE-POINT MEASUREMENTS) 

• SUMMARY: TECHNIQUES UNDER INVESTIGATION 

STREAK RECORDING 
MIE SCATTERING 

LASER-INDUCED PHOSPHORESCENCE 
PULSED GRID RECORDING 

♦DOPPLER MODULATED FLUORESCENCE (DMF) 

• EXPERIMENTAL RESULTS (DMF) 

• RELATED RESULTS FROM AFOSR RESEARCH 
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STREAK RECORDING 


• SEEDED PARTICLES (DROPS) ARE ILLUMINATED BY CW LASER SHEET 

• HIE SCATTERING GIVES STREAKS ON FILM (RETICON) 

• GRADED INTENSITY BEAM (TRANSVERSE TO SHEET) GIVES POTENTIAL 
FOR 3-D RECORDING, IF PARTICLES UNIFORM IN SIZE 



PULSED GRID RECORDING 

• A PULSED TUNABLE LASER IS USED TO ILLUMINATE A GRID IN THE 
FLOWFIELD 

• LASER-INDUCED PHOSPHORESCENCE (OR CHEMILUMINESCENCE) PROVIDES 
LONG-LIVED RADIATION FROM EXCITED FLUID ELEMENTS 

• MULTIPLE EXPOSURE OF GRID PATTERNS ARE RECORDED ON A SINGLE 
FRAME OF RETICON CAMERA (OR FILM) 

9 CANDIDATE MATERIALS: BIACETYL( DROPLET OR VAPOR) 

N0_2/C0 CHEMILUMINESCENCE 
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DOPPLER-MODULATED FLUORESCENCE 




© L4 



FILM RECORD 


RETICON RECORD 













EXTENSION TO SUBSONIC FLOWS 

• nx WAVELENGTH 

• MONITOR VARIATION IN INTENSITY TO INFER VELOCITY 

• RECORD AT HIGH REPETITION RATE TO DETERMINE U(x,y,t) 



SUMMARY 

DOPPLER-MODULATED FLUORESCENCE 

SIGNIFICANCE: MULTIPLE-POINT MEASUREMENTS IN UNSEEDED FLOW 

STATUS: FIRST GENERATION EXPERIMENTS COMPLETED SUCCESSFULLY 

USING SCANNING CONCEPT, 1 BEAM, SUPERSONIC aOW 

LIMITATIONS: LENGTHY RECORDING TIMES 

LIMITED VELOCITY RESOLUTION 5 M/SEC) 

FUTURE: DEVELOP SUBSONIC FLOW TECHNIQUE USING FIXED 

WAVELENGTH, MULTIPLE BEAMS 
DYE LASER TO OPTIMIZE WAVELENGTH 
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MODELING OF FLUCTUATING MASS FLUX 
IN VARIABLE DENSITY FLOWS 


Ronald M. C. So, H. C. Mongia, and M. Nikjooy 
Department of Mechanical and Aerospace Engineering 
Arizona State University 
Tempe, Arizona 85287 


Abstract 

In practical combustion systems, the assumption of chemical equilibrium 
is not always appropriate. This is especially true in regard to predicting 
various critical combustion processes including ignition, flame stability, 
smoke and gaseous emissions. In order to estimate combustor performance 
accurately, it is important to calculate combustor internal profiles of 
temperature, species and velocities to reasonable accuracy. This, in turn, 
requires realistic assumptions concerning the turbulent transports of heat, 
mass and momentum and their interaction with chemical reactions. 

Conventional combustion models assume the scalar transport process can 
be described adequately by an assumed two/three parameter probability density 
function (pdf) and that the transport of heat and mass is essentially identical. 
These models give reasonably well correlated results for simple flows. However, 
their applications to complex flows involving strong recirculation and swirl 
do not appear promising. Since such flows are of common occurrence in any 
combustion systems, straight forward' extension of conventional combustion 
models to combustor flow calculations is not adequate. 

The present approach is formulated to investigate some of the problems 
that are common to both reactive and non-reactive variable density flows. 

Such problems are mass transport in a variable density flow, validity of 
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constant density turbulence models, relative merits of Reynolds versus Favre 
average and calculated versus assumed pdf for the scalar. The approach solves 
for both Reynolds and Favre averaged quantities and calculates the scalar pdf. 
Turbulent models used to close the governing equations are formulated to 
account for complex mixing and variable density effects. In addition, turbulent 
mass diffusivities are not assumed to be in constant proportion to turbulent 
momentum diffusivities. 

The governing equations are solved by a combination of finite-difference 
technique and Monte-Carlo simulation. Some preliminary results on simple 
variable density shear flows are presented. The differences between these 
results and those obtained using conventional models are discussed. 


148 



MOTIVATION 


• UNDERSTANDING OF TURBULENT 
REACTIVE FLOWS 

• MODELING OF TURBULENT REACTIVE 
FLOWS 

• COMBUSTOR FLOW MODELING 

• DEVELOPMENT OF SUITABLE COMBUSTOR 
DESIGN GUIDES 


TURBULENT REACTIVE FLOWS 



FLOW (VELOCITY, SPECIES 





MODELS) 
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MAIN PROBLEM IN COMBUSTION MODELING 


NON-REACTIVE PLOW MODELING 

• CONTINUITY 

• MOMENTUM (TURBULENT MOMENTUM FLUXES, - pn'.n'. 

• ENTHALPY (TURBULENT HEAT FLUXES, L'n' , ETC.) 

• SCALAR (TURBULENT SCALAR FLUXES, 0'ni , ETC.) 

• FLUCTUATING MASS FLUXES (p'u' , ETC.) 

• QUESTION OF REYNOLDS VS. FAVRE AVERAGE 

• VALIDITY OF CONSTANT DENSITY MODEL FOR 
VARIABLE DENSITY FLOWS. 

REACTIVE PLOW MODELING 

• EVALUATION OF MEAN FORMATION RATES 

- R ($, h, p, ) I R ($, h, p, ) 

- FAST/FINITE CHEMISTRY 

• VALIDITY OF CONSTANT DENSITY TURBULENCE MODEL 
FOR REACTIVE FLOWS 


, ETC.) 
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ASSUMPTIONS 


• FAST CHEMISTRY 

(IF IT MIXES, IT REACTS) 

• ONE STEP FORWARD REACTION 

(FUEL + OXIDANT - PRODUCTS) 

• TURBULENT TRANSPORT OF DIFFERENT 
SPECIES IS THE SAME 

(TURBULENT SCHMIDT NUMBER 
FOR ALL SPECIES IS IDENTICAL) 

• EDDY DIFFUSIVITIES FOR MASS AND 
HEAT ARE IDENTICAL 

(LEWIS NUMBER EQUALS UNITY) 


CONSERVED SCALAR, f 

f = MIXTURE FRACTION 

= MASS OF SPECIES 1 
MIXTURE MASS 

RELATED TO OTHER CONSERVED SCALAR, 
SUCH AS ENTHALPY, h, BY 

^ ^ox 
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EQUATIONS SOLVED IN CONVENTIONAL REACTIVE 
FLOW MODELING 

. MEAN CONTINUITY 

• MEAN MOMENTUM (CLOSURE VIA TURBULENCE 
MODELING) 

. MEAN CONSERVED SCALAR (CLOSURE VIA 
TURBULENCE MODELING AND SCHMIDT NO.) 

. MEAN VARIANCE OF CONSERVED SCALAR 
(f^, WHERE f' = f-f) 


CONVENTIONAL APPROACH 


• SOLVE EQUATIONS USING FINITE 
DIFFERENCE SCHEME 


• ASSUME A TWO-PARAMETER PDF 
FOR f TO ACCOUNT FOR SCALAR 
FLUCTUATIONS 


• DETERMINE THE PARAMETERS BY 
REQUIRING FIRST AND SECOND 
MOMENTS TO AGREE WITH f AND 
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PRESENT MEASUREMENT AND PREDICTION CAPABILITIES FOR TURBULENT REACTIVE FLOWS 




Measurement Capabilities 

Prediction Capabilities 



Physical 

Quantities 

Probe 

Data 

Optical 

Data 

Difficult 

? 

Assumed 

Pdf 

Pdf 

Transport 

Difficult 

7 

Remarks 

Flow 

Field 


X 



X 





X 

X 


X 



Hot-wlre/LDA 

technique 



X 

m 

X 

v.w u 




Hot-wlre/LDA 
technlque;Need 
to Improve model 


Hi 




m 


X 

n 


h'u; 

X 



X 


X 

ti 


Pdf of 

X 


X 

(Major 

Species) 

X 

X 


Raman scattering 
tech, need to 
Improve model 

Combustion 

Field 



■ 


X 

(Model ) 

X 

(Model) 

X 



Bi 


X 

X 

X 

1 

1 

Deduced from 
meas,; Rayleigh 
scattering 


fm 


X 

(N2) 

X 


X 


Raman Scatter- 
ing. O 9 diffi- 
cult 

Pollutants 

0A T'. Etc 

U2 


X 

(N2) 

X 


X 


n 





X 
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CONVENTIONAL APPROACH TO 
COMBUSTION MODELING 


• REDUCE PROBLEM TO EQUIVALENT 
SCALAR MIXING PROBLEM 











































LIMITATIONS OF CONVENTIONAL APPROACH 


. INFINITELY FAST CHEMISTRY 

. ONE/TWO STEP REACTIONS ONLY 

. SCALAR FLUCTUATIONS DETERMINED BY 
ASSUMED PDF 

. ALL SPECIES DIFFUSE AT SAME RATE 
• HEAT/MASS DIFFUSE AT SAME RATE 
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PRESENT APPROACH ADDRESSES THE 


FOLLOWING PROBLEMS 

• PDF OF SCALAR 

• MASS DIFFUSION IN TURBULENT FLOWS 

• REYNOLDS OR FAVRE AVERAGE 

• CALCULATION OF FLUCTUATING 
MASS FLUXES 


NOMENCLATDRE 

, etc. 

etc. 

, 0, etc. 

* 0*» P’ * etc. 

F (0 ) 
a 

< > 


FAVRE AVERAGE 
REYNOLDS AVERAGE 

FLDCTOATING PART OF FAVRE DECOMPOSITION 

FLDCTDATING PART OF REYNOLDS 
DECOMPOSITION 

PDF. OF 0 
TIME AVERAGE 
INSTANTANEOUS QUANTITIES 
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GOVERNING EQUATIONS 


• CONTINUIIT 


d 

/ — « 


P 

dx. 

\ 

i 


MOMENTDM 

d 

r " „ 


p 

dx . 

L 

1 


PDF OF 

Mai 


3P 

dx 


U. 

1 



^13 

rj,aF(?„) • 

1 - — f 

p dx. 
^ 1 

L dx. 
1 

1 d0 L 

a 


F(0 ) S(0 ) 


+ E (0 , X.) 

a 1 
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TDRBDLENT CLOSURE 


TWO-EQUATION PLUS NOtJ-EQUILIBRlUM ALGEBRAIC STRESS MODEL 


P U, 


3k a fr k _ , . ak 1 

= |C, _ p <u. u.> I 

ax. dx. L e ^ dx. J 


- p <u u > iHi 

^ ax. 


- pc 


- ^ <P^i> ap_ 

p ax. 


— tt 3 To ✓ V 3c 1 

P U^ = Cg _ p <u n > - 

aXj 3x^ L e ^ dx, J 


e — 


k ^ J 0Xj k 


. e <p'ni> aP ^ _ aUi 

+ !- — I + pe 


k p ax. 


ax. 


or- 


[<"‘ ^ i? ] 




- <U, U_> 1 + < £ / + 


i m 


££i 1 

ax J 

m 


3 ua 9 u 


\— =■ ^ — <-) > 

p 'ax. dx. 

j 1 
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MODEL FOR p AND <p'u^ 


CONSIDER MIXING OF TWO FLUIDS WITH DIFFERENT DENSITIES. 


P 

U. = U. + a.^<u.0> 

1 1 2 1 


<p'u*> = - a. p <u.0> 

1 2 1 


1 



$ FROM PDF F(0 ) AND <n.0> FROM TURBULENCE MODEL 


MODEL FOR <u^0> 


NON-EQUILIBRIUM ALGEBRAIC FLUX MODEL 


3$ dU^ p d0 

- <n^ Uj> - <Hj0> + < - > r 0 


dz . 
J 


dz 


j 


P dz. 
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SOLUTION PROCEDURE 


• FINITE DIFFERENCE SOLUTION OF ALL GOVERNING EQUATIONS 


• OPERATOR-SPLITTING TECHNIQUE APPLIED TO PDF EQUATION 


ASSUME F(0 ) = F*(P ) + F'(0 ) 

a a a 

F* AND F' GIVEN BY SOLVING 


6F* ^ 1 r £F* 1 

^ 3x. p dx, L dx J 

1 i i 


dF* 

dt 


E = 


2“ \ { F* 

- 2Pj 



(0'^) p* (20 

F* (0 ) 
a 


+ E(0 , X.) 

a 1 


- 0'^) d0' 


a 
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MONTE CARLO MODELING IN ELLIPTIC FLOWS 


Krishnan Radhakrishnan 

Department of Mechanical Engineering and Applied Mechanics 
The University of Michigan 
Ann Arbor, Michigan 48104 


A Monte Carlo method capable of predicting scalar (e.g. concentration) prob- 
ability density functions (pdf’s) in nonreacting and reacting elliptic flows 
is being developed. In particular, attention is being focused on a research 
combustor specifically designed to provide measurements of concentration and 
temperature pdf’s in nonreacting and reacting flows. (e.g. 8). 

The modeling technique being developed involves the coupling of a Monte 
Carlo coalescence/dispersion model (e.g. 2-7) with currently available codes 
(such as the finite-difference scheme described in (1)) for elliptic flow- 
field calculations. The latter provides solutions for the mean velocity, 
turbulent kinetic energy and dissipation rate. These serve as Inputs to 
the Monte Carlo calculation procedure which is used to solve the transport 
equations for scalar pdf’s. 

The attraction of the Monte Carlo method is its practicability for 
multi-dimensional pdf’s, while conventional finite-difference solutions of 
the pdf transport equations require prohibitive amounts of computer time 
(e.g. 2). This method also eliminates the necessity of assuming a shape 
for the pdf; instead, given an initial distribution function, the model 
allows for the temporal evolution of the pdf. Also, it can account for 
finite rate chemical reactions, which is essential for the successful 
prediction of pollutant formation rates and flame stability. 

In the Monte Carlo calculation procedure, the composition distribution 
(due to concentration fluctuations) at a point in the flowfield is described 
by a statistical ensemble of computational elements which have a distribution 
in thermodynamic states. Each element has scalar properties associated 
with it so that its thermodynamic state is well defined. The number of such 
elements in a given state can vary because of turbulent mixing, chemical 
reactions, etc. Given such an ensemble of elements at every grid point in 
the flowfield, mean values of any scalar, as well as variance or higher 
moments can be obtained from instantaneous averaging over the ensemble pop^ 
ulation. 

The effects of mixing, convection, diffusion, and chemical reaction on 
the scalar pdf’s are incorporated into the calculation procedure as follows. 

The turbulent and moecular mixing processes are represented by allowing 
elements within each ensemble to mix with one another. These mixing inter- 
actions occur at discrete time intervals given by a mixing frequency which 
can be related to the local turbulent kinetic energy and dissipation rate 
(see reference 2 for details) . The effects of turbulence on mixing and 
chemical reaction are thus introduced through the mixing frequency. During 
the time between mixing interactions, each element undergoes chemical 
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reaction at a rate specified by the user. Either an infinite rate assumption 
can be made or a chemical reaction mechanism and rate data specified^ Since 
the processes of mixing and reaction are treated separately, the use of 
multiple finite rate reactions does, not increase the complexity of the form- 
ulation. Also, the instantaneous reaction rate for each element is calculated, 
so the mean reaction rate can be obtained by simply averaging over the ensem- 
ble population. The transport processes of convection and diffusion are also 
treated separately; the former by transporting elements from one ensemble 
(i.e. grid point) to another at a rate governed by the local fluid velocity, 
and the latter by exchanging elements between ensembles at adjacent nodes at 
a rate governed by the local turbulent diffusivity. 

The specification of initial ensembles at the grid points is as follows. 
For elliptic flows, properties at a point can be affected by values down- 
stream. The calculation procedure will therefore involve iteration on some 
assumed initial distributions. The processes of convection, diffusion, 
mixing, and reaction will be simulated as described above. Also, the inlet 
flow conditions will specify ensembles along one axial location. The itera- 
tion will continue until the ensembles have become time stationary; a 
correct ensemble will then have been created at all spatial locations in 
the flowfield. 

The ability of the Monte Carlo method to provide the coupling of finite 
rate fluid mechanic mixing and finite rate combustion chemistry is illustrated 
by presenting a coalescence/dispersion model developed for studies of ignition 
and blowout in a combustor primary zone. The use of the flame stability model 
(described in references 5-7) is demonstrated by generating lean ignition 
limits of premixed turbulent flames. The variations of the predicted lean 
ignition limit with reference velocity, inlet temperature, and pressure 
compare favorably with experimental trends and magnitudes. 


This work is supported by NASA Grant NAG 3-294. 
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Method 


Combine existing computer codes for elliptic flow field calcu- 
lations with a Monte Carlo method to simulate the transport 
equations for scalar probability density functions. 
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Solution Technique For Elliptic Flows 


Interative as follows: 


Assume initial ensemble (and hence, pdf's) at all grid points. 
Calculate < p >, < T >, etc. 

Solve for flow velocity, turbulent kinetic energy and dissipa- 
tion rates at all grid points. 

Generate new pdf's at all grid points. 

Check if new pdf's agree with old values. 


Schematic of Bluff Body Slobilized Flome 


Recirculotion Zone 
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Simulation of Recirculation Zone 


Unburnt 
Mixture ^ 


Porfiolly Stirred Reoctor 



Average 
M ixture 


Monte Carlo Mixing Interaction 



• 

f 

o 


E lement 

i 

j 

‘ j 

Burnt 

Froction 

Bi 


t Bj + B. ) / 2 
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Model Parameters and Relations 


Mixing Frequency 

e = Cg (pyML^)^/^ 

Inflow Time 

T^. = (V/m.v) 

Residence Time 

X = (V/ifigV) 

Chemical Kinetics 

Two Step Mechanism 

C^Hg + 3.5 O 2 3 CO + 4 H 2 O 
CO + 0.5 O 2 t CO 2 

Reaction Rates from Westbrook and Dryer (1981) 
Preexponential Constants Adjusted 

Ensemble Property 


166 



Ensemble Burnt Fraction 


Example : Ignition 


f(B) 



Development of 
Distribution Function 
Due to Mixing 


f(B) 



Ignition 
B Increases 



B 


No Ignition 
B Decreases 


EXAMPLE .LEAN IGNITION LIMIT 



( t/x:) 


O 












Variation of Predicted Limit with Pressure 
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COMPUTATIONS OF EMISSIONS USING 
A 3-D COMBUSTOR PROGRAM 


S. K. ’Srivatsa 

Garrett Turbine Engine Company 
Phoenix, Arizona 85010 

Objective; Extend a 3-D combustor program; 

• To predict pollutant emissions of soot and NO^ 

• To include the influence of soot, CO 2 , and H 2 O on radiation heat 
transfer 

• To extend the two-step hydrocarbon oxidation mechanism to a more 
detailed four-step scheme 

Background ; EPA regulations for aircraft gas turbines make it imperative that 
soot and NO emissions be controlled. The use of alternative synthetic 

A 

fuels, recently being considered, results in significant increases in 
soot production. Thus, an improved understanding of the physical and 
chemical processes governing soot production is necessary. Since these 
processes are not well understood, only an approximate modeling of soot 
production is currently possible. Lack of general 3-D computer codes has 
also hindered research in this area. In the present program, soot and 
NO emissions were computed by using an existing general 3-D computer 
program. This will significantly aid in the development of soot and 
NO production models in practical systems. 

Approach; A general 3-D combustor performance program developed by Garrett 
was extended to predict soot and NO^ emissions. The soot formation and 
oxidation rates were computed by quasi-global models, taking into account 
the influence of turbulence. Radiation heat transfer was computed by the 
six-flux radiation mode. The radiation properties include the influence 
of CO 2 and H 2 O in addition to soot. NO^ emissions were computed 
from>a global four-step hydrocarbon oxidation scheme and a set of 
rate-controlled reactions involving radicals and nitrogen oxides. 
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Results: Computations performed for a plug flow reactor show the four-step 

scheme to be far superior to the two-step scheme in predicting 
temperature and species concentrations. 

• Computations were performed for idle, cruise, and takeoff conditions of 
a JT8D combustor. These showed that the present model is capable of 
producing reasonable predictions of smoke and NO emissions and of 

A 

the wall radiation flux. 
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COMPUTATIONS OF EMISSIONS USING 
A 3-D COMBUSTOR PROGRAM 


OBJECTIVE: 

EXTEND AN EXISTING 3 -D COMBUSTOR PROGRAM; 

• TO PREDICT POLLUTANT EMISSIONS OF SMOKE AND NOx: 

• TO INCLUDE THE INFLUENCE OF SOOT, CO2. AND H2O ON 
RADIATION HEAT TRANSFER; AND 

• TO EXTEND THE TWO-STEP HYDROCARBON OXIDATION 
MECHANISM TO A MORE DETAILED FOUR-STEP SCHEME 


APPROACH 

PROGRAM INVOLVED FOUR TASKS 

• TASK I - FORMULATION OF THE METHOD 

• TASK II - COMPUTER CODING 

• TASK III - COMPUTATION OF TEST CASES 

> IDLE. CRUISE. AND TAKE-OFF CONDITIONS FOR A 
JTBD COMBUSTOR 

• TASK IV - REPORTING AND DOCUMENTATION 
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EXISTING 3-D COMBUSTOR PROGRAM 


• GENERAL PROGRAM FOR 

• RECIRCULATING. SWIRLING. TURBULENT. REACTING FLOW 

• GASEOUS AND/OR LIQUID FUEL COMBUSTION 

■ RADIATION HEAT TRANSFER 

• GENERAL TRANSPORT EQUATION 

div b>U^ • GRAD <(>] « S 0 

• VARIABLES: ITp. k. r. 1 . M|u. Mco (Mq,. Mco?. MhoQI 
It IT. p]. Rx. Ry. Rj SPRAY DYNAMICS/COMBUCTION 

• PHYSICAL MODELS: 

• TURBULENCE: k-c MODEL 

• CHEMISTRY: 2 -STEP REACTION SCHEME 

• CHEMICAL REACTION RATE: MODIFIED EBU MODEL 

• RADIATION: SIX-FLUX MODEL 

• NUMERICAL 

■ FINITE-DIFFERENCE ITERATIVE METHOD 

■ SUITABLE FOR COMPLEX GEOMETRIES 

• SUITABLE FOR NON-UNIFORM GRID SPACING 


FOUR-STEP HYDROCARBON OXIDATION MECHANISIVI 


• TWO-STEP SCHEME IN ORIGINAL 3 D PROGRAM: 

CxHy + O2 - CO + H2O 

CO + D2 “* CO2 

• FOUR-STEP SCHEME: 

CxHy — CxHy.2 + H 2 
CxHw -2 + O2 - CO + H 2 
CO + O2 - CO2 

H2 + O2 — H2O 

• FOUR-STEP SCHEME DESCRIBES THE FOLLOWING BASIC STEPS 
OF HYDROCARBON OXIDATION 

• TRANSFORMATION OF HYDROCARBON FUEL INTO INTERMEDIATE 
HYDROCARBONS AND HYDROGEN WITH LITTLE RELEASE OF ENERGY 

> OXIDATION OF INTERMEDIATES TO CO AND H2 

> OXIDATION OF CO TO CO2 

> OXIDATION OF H2 TO H2O 

• TWO ADDITIONAL EQUATIONS FOR: CxHy .2 AND H2 
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NOx EMISSIONS 

• SPECIES CONSIDERED: CxHy. CxHy. 2 , CO. N. NO. NO2. H. 

H2. oh. 0. CO 2 . H 2 O. O 2 . N2 

• TRANSPORT EQUATIONS FOR: CxHy. CO. CxHy. 2 . N. NO. NO2. 
H. H 2 . oh. 0 

• CO2. H2O. O2 FROM C. H. 0 ELEMENT CONSERVATION 

• Mh2 = 1-lM 

• SOURCES OF SPECIES COMPUTED FROM MODIFIED EBU MODEL 

• CxHy. CxHy. 2 . H 2 . AND CO SOURCES FROM 4 -STEP 
REACTION SCHEME 


NOx emissions 


• SOURCES OF N. NO. NO2. H. H2. OH. 0 FROM REACTION MECHANISM: 


0 + N 2 = N + NO 
NO 0 = O2 + N 
N + OH = NO + H 
N 2 + O2 = N + NO2 
N 2 + O2 = NO NO 
NO + NO = N -I- NO2 
M■^N 0 = M + N + 0 
M + NO2 = M + 0 + NO 
M + NO2 = O2 N 4 M 
NO 4 O2 = NO2 4 0 


H 4 NO2 = NO 4 OH 
H2 O2 “ OH 4 OH 
OH 4 H2 = H2O 4 H 
OH 4 OH = 0 4 H2O 
0 4 H 2 = H 4 OH 
H 4 O2 = 0 4 OH 
M404H = 0H4M 
M 4040 = 024 M 
M 4 H 4 H=H 24 M 
M 4 H 4 OH = H2O 4 M 


STRONGLY COUPLED NONLINEAR STIFF EQUATIONS 

PRATT'S CREK PROGRAM FOR STIFF EQUATIONS ADAPTED TO 
3-0 COMBUSTOR PROGRAM 
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SOOT EMISSIONS 


• VARIOUS STEPS INVOLVED IN SOOT FORMATION/OXIDATION NOT QUANTIFIED 

• HENCE QUASI GLOBAL MODELS CONSISTING OF THREE STEPS 

■ SOOT miEI FORMATION 

> SOOT PARTICLE FORMATION 

■ SOOT PARTICLE OXIDATION 

• TRANSPORT EQUATIONS FOR NUCLEI AND SOOT CONCENTRATIONS WITH 
SOURCES CONTAINING REACTION RATES 

• REACTION RATES CONTAIN FORMATION AND OXIDATION RATES DEPENDENT ON: 

■ UNBURNT FUEL AND O2 CONCENTRATION 

■ TEMPERATURE 

• TURBULENCE PARAMETERS c/k (SIMILAR TO EODY-BREAK-UP MODEL) 

• TWO SIZES OF SOOT PARTICLES CONSIDERED: ONE EQUATION FOR EACH SIZE 


RADIATION HEAT TRANSFER 

• SIX-FLUX RADIATION MODEL - SCHUSTER AND HAMAKER 

• MAJOR CONTRIBUTORS TO RADIATION = SOOT. CO2. H2O 

• RADIATION PROPERTIES DEPENDENT ON SOOT. CO2. H2O CONCENTRATIONS 
f TOTAL PROPERTIES (INTEGRATED OVER WAVELENGTH| 

• CO2 - H2O EMISSIVITY CORRELATIONS INCLUDING OVERLAP AND 
PRESSURr CORRECTION FACTORS 

eC+W = fC + 

»C+W = Ln( 1 - iQ+vvl ! S 

• SOOT RAOIATION PROPERTIES 

■ ABSORPTION COEFFICIENT COMPUTED FOR BOTH PARTICLE SIZES 
AND ADDED TOGETHER 

• GAS-SOOT MIXTURE. 

■ ABSORPTION COEFFICIENTS OF GAS AND SOOT ADDED 
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RESULTS 


• FOUR-STEP SCHEME 

■ PLUG FLOW REACTOR FOR LEAN. STOICHIOMETRIC. ANO 
RICH PROPANE FLAMES 

■ CONSIDERABLE DEVIATIONS OF TWO-STEP PREDICTIONS 
FROM MEASUREMENTS 

■ FOUR-STEP PREDICTIONS AGREE CLOSELY WITH 
MEASUREMENTS 

■ SOME DISCREPANCY IN FOUR-STEP H2 PREDICTIONS AT 
HIGHER EQUIVALENCE RATIOS 


LEAN C 3 HB FLAME (0 = 0.12) 
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TEMPERATURE. K 


LEAN C 3 H 8 FLAME [0 




AXIAL DISTANCE. CM 


AXIAL DISTANCE. CM 


LEAN C 3 H 8 FLAME (0 = 0.12) 


\\ 

\c,H. WTIRMfOIATtf \N 


AXIAL DISTANCE. CM 



AXIAL DISTANCE. CM 







JT8D-17 COMBUSTOR 












PREDICTED EMISSIDNS INDEX WITH FDUR-STEP 
HYDRDCARBDN DXIDATIDN SCHEME 



EMISSION INDEX 
Gm OF EMISSIONS/Kg OF FUEL 


CONDITION 

SMOKE 

NOx 

IDLE 

0.056 (0.6) 

'=0 

CRUISE 

1.3 

13 

TAKEOFF 

]2 |2.8) 

27 (24.4) 


NOTE: VALUES IN PARENTHESES ARE EXPERIMENTAL MEASUREMENTS 

• PREDICTED RADIATION FLUXES IN THE SAME RANGE AS 

MEASUREMENTS FOR SLIGHTLY DIFFERENT OPERATING CONDITIONS 


CDNCLUSIDNS 


• MODEL PRODUCES REASONABLE RESULTS FOR EMISSIONS 
AND RADIATION FLUX 


• LACK OF EXPERIMENTAL DATA PRECLUDES MORE DETAILED 
MODEL VALIDATION ' 


181 




FAST COMPUTATIONAL KINETICS PROGRAM 


David T. Pratt 

University of Washington FU-lO 
Seattle, Washington 98195 

1 . Introduction 

This project has two goals: First, to develop very fast algorithms 

for solution of the transient and stationary-state equations describing 
combustion chemical kinetics in batch and stirred-tank chemical reactors, 
for imbedding in larger hydrodynamic codes which model mixing, heat 
transfer and chemical change in gas turbine combustors and other combus- 
tion devices. The second goal is to develop "user-friendly" interactive 
computer programs which combine the fast kinetics algorithms with graphics 
and user selection menus to create a dynamic environment for parametric 
variation, preliminary design, and in general permitting the user to gain 
an improved understanding of the thermal -chemical -phys ical behavior of 
combustion processes, without being hampered by considerations of numerical 
methods, stability, convergence or computational efficiency. 

2. Batch Kinetics (l-D) Algorithm Development 

A typical batch combustion problem at constant pressure consists of 
three distinctly different chemical-physical regimes, each of which 
requires a unique algorithm for the peculiar form of stiffness of the 
governing set of ordinary differential equations. (The term "stiff", 
coined by Hi rschfelder , means that the characteristic time constants of 
the coupled set of equations differ by many orders of magnitude, so that 
integration of the equations by conventional methods would require time- 
steps very much smaller than desired, resulting in excessively large 
computational work. [1]) 

A single-step integration algorithm was developed which automatically 
identifies equations as having either positive or negative time constants, 
and then approximates the unstable equations (positive time constants) 
with the trapezoidal rule, and the stable equations (negative time 
constants) with a decaying exponential function. This "exponent ial -fi tted 
trapezoidal rule" was adapted form the work of Liniger and Willoughby [2] 
and of Brandon [3] • 

During induction and early heat-release regime, the species equations 
are dominated by positive time constants, and the temperature also 
exhibits a positive time constant. Since very smal 1 steps are required 
for integration of unstable equations, a simple predictor-corrector scheme 
with functional iteration assures the least computational work possible. 
However, during late heat release and equilibration, when the temperature 
and species equations exhibit negative time constants, large stepsizes 
are now possible, so that Newton-Raphson iteration with calculation of 
the full Jacobian matrix is the optimal convergence method. At the present 
level of development, the batch kinetics code CREK-ID (for "combustion 
reaction kinetics, one-dimensional") appears to be at least five times 
faster than LSODE, the Lawrence Livermore Laboratory implementation of 
the Gear-B algorithm [4]. Further details of the algorithm are given 
in Reference 5. 
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3. Stirred Reactor (O-D) Algorithm Development 


An earlier packaged code, CREK [6], solved for time-stationary 
states of a perfectly-stirred chemical reactor. A major improvement 
has been made in the automatic estimate-generating feature of the CREK 
code. In this approach, the transient sti rred-reactor equations are 
solved for the predictor step only, using the same techniques developed 
for the batch kinetics CREK- ID code. The resulting approximate solution 
gives a very good estimate of the steady-state solution, which is then 
solved for exactly using the original CREK algorithm. As a result, the 
speed and reliability of the "CREK-OD" program appears to have been 
increased by at least an order of magnitude over the original CREK code. 

4. Interactive Computer Programs 

Two "user-friendly”, menu-driven programs have been written to 
date. The first program, EQLBRM, solves for chemical equilibrium states 
for combustion of arbitrary hydrocarbon fuels with air, at a specified 
constant pressure, fuel-air ratio and initial temperature of fuel and 
air. The second program, AVC0-M2, is a reactor-theoretic simulation of 
a gas turbine combustor [7]. The combustor is assumed to consist of 
up to nine "flow elements" connected in series, and optionally one flow 
element connected in recycle. Each flow element may be a zero-dimensional 
or perfectly-stirred reactor (PSR) , a one-dimensional or plug-flow 
reactor (PFR) , or a non-reacting mixer (MIX). The user is prompted to 
select the order, type and size of the flow elements, the recycle path, 
and the distribution of fuel and air into each individual flow element. 

On demand, the program also displays a crude graphic schematic of the 
combustor model under consideration. 
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Variation with time of temperature and chemical species 
mole fraction in adiabatic batch reaction. 
(Pyrolized CH^-air, Stoichiometric, T=1000K, P=10 atm) 
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Conventional staggered finite-difference grid for calculation 
within a flow domain. Points H (“high") and L ('-'low") are 
above and below the plane of E-W-N-S-P. PP refers to condi- 
tions at P at the previous time step. 
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HOMOGENEOUS, ADIABATIC, GAS-PHASE 
REACTION AT CONSTANT PRESSURE 


1 



_L 4. JL rr ^ 

■ P PV 


Batch kinetics (1-D) 
Stirred-reactor kinetics (O-D) 

where 


r, (a:j - crjHRj - R.j) 


^ N 

R. = k. (paj J n (pa. ) 

j J m k=l " 


kj 


k. = lO^j T^j exp(-T./T) 

J J 


p = P/RTa, 


m 


N 

= E a- 
k=l ' 


N 

H„ = Z h.a. = const. 

0 k=l k !< 


3 

1 
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"LOCALLY EXACT" SOLUTIONS 
OF THE TRANSIENT EQUATIONS 


1 . Decoupl i ng 




A. . B.a. 


2. Integration 


n+1 


“i 




a" . 

“1 ^ B7 




"+1 

or a, ' Of + hfj 


B,h 


where h = t'^''’^ - t^ 


3. Analysis 

The locally exact solution is stable (for h — > “) if and 
only if B^.< 0. 
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FORMALLY DECOUPLING THE 
EOUATIONS: "L - FORMULATION" . 


let 






- Q. 




then 


dt 


q, - L,a. . i (o, - a,) 


Q. tio* 
pV 1 


I + 

“i 


''i = 


f 


ni * 

Q • “*■ “TT O . 

pV 1 






. . .A^ is positive-definite; negative-definite 
*. "L-formulated" equations are always 
stable , but may be inaccurate. 
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FUNCTIONALLY DECOUPLING THE 
EQUATIONS : "Z-FORMULATION" 


let 


da. 


= f.. = 


dt 'i 


xO . 7O / 0\ 

f i * Z, (o, - a, ) 


where Z. e 


df. 


i da^ 


Z + -?r 


k=1 H 




(and where t = - Z h|^f|^ / Z C„ a^) 


k Pk 


da. 


then 


' - [f? - Z°a°] + Z°a. 


dt 


=> A. = 


[f“ 


7O On 

Z,o,] 


6 . 

1 




. . . neither A^. nor are positive-definite or negative-definite. 

"Z-formulated" equations are always accurate , but may be unstable . 
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HARK2I )K>Jo»c 

■MARK2I* IS AN INTERACTIVE VERSION OF THE MARK-II COMBUSTOR MODEL* 
THIS IS A PRELIMINARY DESIGN TOOL, A MEANS OF GAINING INTUITIVE INSIGHT 
INTO EFFECTS OF CHANGES IN FUEL-AIR MIXING OR PARTITIONING ON TURN-DOWN 
RATIO, COMBUSTION EFFICIENCY AND POLLUTANT FORMATION RATES* 

AN INITIAL DATA SET IS TAKEN FROM DATA FILE -MARK2*DAT- BUT CAN BE ALTERED 
INTERACTIVELY, AND USED IN CONSECUTIVE RUNS* 

MAF^K-II REPRESENTS A SIMPLE BRAGG COMBUSTOR CONSISTING OF A 
MAXIMUM OF 9 FLOW ELEMENTS WITH THE ADDITION OF A SINGLE RECYCLE 
ELEMENT* FLOW ELEMENT TYPES MAY INCLUDE: 

1) NON-REACTING MIXERS CMIX*), IN WHICH THE CHEMICAL REACTIONS ARE 
ASSUMED TO HAVE STOPPED DURING THE MIXING PROCESS; 

2) PERFECTLY STIRRED REACTORS CPSR*), WITHIN WHICH INTENSE SELF- OR 
DACK-MIXING IS ASSUMED TO OCCUR, SO THAT THERE ARE NO AXIAL GRADIENTS; 

3) PLUG FLOW REACTORS CPFR*)* 

THE USER MAY DEFINE THE MODEL AS HAVING UP TO 9 ELEMENTS IN 
SERIES WITH AIR AND FUEL INLET JETS AT EACH ELEMENT* THE RECYCLE ELEMENT 
MAY DE OF ANY OF THE THREE FLOW TYPES, AND MUST RECYCLE FROM A HIGHER 
NUMBERED ELEMENT TO A LOWER* COOLING BOUNDARY LAYER EFFECTS AND CHEMICAL 
REACTIONS WITHIN THE BOUNDARY LAYER ARE NOT CONSIDERED* 

PLEASE WAIT A MOMENT WHILE INITIALIZATION IS COMPLETED. 

— INITIALIZED ~ PRESS <RETURN> TO BEGIN — 


*** INPUT DATA 


FLOW 

AREA 

LENGTH 

FLOW 

ELEMENT* 

(SQ*IN) 

(INCHES) 

TYPE 

1 

1*4A00ET02 

4*0000E-01 

PSR 

n 

l*4600E+02 

1*0000E-01 

MIX 

3 

1*4A00ET02 

l*5000E+00 

PSR 

4 

1*4S00EF00 

2*0000E-01 

MIX 


INLET AIR 
<LBM/S) 
l*3450E+00 
6*9600E-"01 
0*0000E-01 
7*9500E-01 


INLET FUEL 
(LBM/S) 
9.7200E-0 
O.OOOOE-0 
0*0000E-01 
0*OOOOE-Ol 


AIR TEMP « 2»1000ET02 F 

FUEL TEMP 8*0000E+01 F 


COMBUSTOR PRESSURE = 2*2100E+00 ATM 

LOWER HEATING VALUE ^ 1*8500E+04 BTU/LBM 


SELECT AN OPTION BY NUMBER: 

-0- RUN WITH THIS DATA SET 
-1- CHANGE AIR TEMPERATURE 
-2“ CHANGE FUEL TEMPERATURE 
•3- CHANGE LOWER HEATING VALUE 
OPTION? (0-7) 7 


-4- CHANGE NOMINAL COMBUSTOR PRESSURE 
-5- CHANGE RECYCLE ELEMENT STATUS 
-6- CHANGE FLOW ELEMENTS STATUS 
-7- INSPECT SCHEMATIC MODEL LAYOUT 


MARK-II MODEL SCHEMATIC LAYOUT 


A F A A 

1 ! / 1 ! 

^ 

* 1 t t 2 t t Z t ^ A t 
->*PSR)|t— >*MIX* — >»PSR3|c— »»cMIX5»c — > 


PRESS <RETURN> TO CONTINUE 
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%%% INPUT DATA 


FLOW 

AREA 


LENGTH 

FLOW 

INLET AIR 

INLET FUEL 

ELEMENT* 

(SQ*IN) 


(INCHES) 

TYPE 

(LBM/5) 

(LBM/S) 

1 

1 ♦>>600ET02 


>^.0000E-01 

PSR 

l*3A50E+00 

9t7200E-02 

n 

1.4600E+02 


1*OOOOE-Ol 

NIX 

6*?600E-"01 

O.OOOOE -o:i 

3 

1.4600E+02 


l»5000E+00 

PSR 

0,0000E“01 

0*0000E-01 

A 

l*4500E+00 


2*0000E‘ 01 

MIX 

7.9S00E--01 

O.OOOOE-01 

AIR TEMP = 

2*1000ET02 

F 

COMBUSTOR 

PRESSURE 

^ 242:l00E+00 

ATM 

UEL TEMP 

OiOOOOE+Ol 

F 

LOWER HEATING VALUE 

= l*8500E+0^ 

BTU/LBM 


SELECT AN OPTION BY NUMBER; 

-0- RUN WITH THIS BATA SET 
-1- CHANGE AIR TEMPERATURE 
-2- CHANGE FUEL TEMPERATURE 
-3- CHANGE LOWER HEATING VALUE 
OPTION? <0-7> 5 


-A- CHANGE NOMINAL COMBUSTOR PRESSURE 
-5- CHANGE RECYCLE ELEMENT STATUS 
-A- CHANGE FLOW ELEMENTS STATUS 
-7- INSPECT SCHEMATIC MODEL LAYOUT 


RECYCLE ELEMENT STATUS: 

CROSS-SECTIONAL AREA - I.4600EL02 SO. INCH 

ELEMENT LENGTH ^ l.OOOOE+00 INCHES 

FLOW TYPE MIX 

RECYCLES 0.00% OF ELEMENT * 1 OUTFLOW TO ELEMENT ♦ 1 INFLOW 

SELECT AN OPTION BY NUMBER; 

0- NO MORE CHANGES TO RECYCLE ELEMENT 

1- CHANGE AREA 

2- CHANGE LENGTH 

3- CHANGE FLOW TYPE 

A- CHANGE RECYCLE INTAKE PERCENTAGE 
5- CHANGE RECYCLE PATH 
ELIMINATE RECYCLE 

OPTION? <0-6) A 

ENTER PERCENT OF OUTFLOW TO BE RECYCLED: 20, 


RECYCLE ELEMENT STATUS; 

CROSS-SECTIONAL AREA “ 1.4600E+02 SO. INCH 

ELEMENT LENGTH - l.OOOOE+00 INCHES 

FLOW TYPE - MIX 

RECYCLES 20.00% OF ELEMENT * 1 OUTFLOW TO EL.EMENT ♦ 1 INFLOW 

SELECT AN OPTION BY NUMBER; 

0- NO MORE CHANGES TO RECYCLE ELEMENT 

1- CHANGE AREA 

2- CHANGE LENGTH 

3- CHANGE FLOW TYPE 

4- CHANGE RECYCLE INTAKE PERCENTAGE 

5- CHANGE RECYCLE PATH 

6- ELIMINATE RECYCLE 

OPTION? (0-6) 5 

ENTER ELEMENT WHOSE OUTFLOW IS TO BE RECYCLED: (1- 4) 3 
ENTER ELEMENT WHOSE INFLOW IS TO RECEIVE RECYCLE; (1-3) 2 


191 


RECYCLE ELEMENT STATUS: 

CROSS-SECTIONAL AREA 1.^600E+02 SQ/INCH 
ELEMENT LENGTH = l,OOO0E+OO INCHES 

FLOW TYRE == MIX 

RECYCLES 20*00% OF ELEMENT ♦ 3 OUTFLOW TO ELEMENT ♦ 2 INFLOW 

SELECT AN OPTION BY NUMBER: 

0- NO MORE CHANGES TO RECYCLE ELEMENT 

1- CHANGE AREA 

2- CHANGE LENGTH 

3- CHANGE FLOW TYPE 

4- CHANGE RECYCLE INTAKE PERCENTAGE 

5- CHANGE RECYCLE PATH 

6- ELIMINATE RECYCLE 

OPTION? (0-6) 0 


INPUT 

DATA *** 






FLOW 

AREA 

LENGTH 

FLOW 


INLET AIR 

INLET FUEL 

ELEMENTI 

(SQ*IN) 

(INCHES) 

TYPE 


(LBM/S) 

(LBM/S) 

1 

l*4600E+02 

4.0000E-01 

P9R 


l*3450Ef00 

9*7200E-02 

2 

1*4600ET02 

l.OOOOE-01 

MIX 


6*9600E-01 

0*0000E-01 

3 

l*4600E+02 

1 ^SOOOETOO 

PSR 


0*OOOOE -01 

0*0000E-01 

4 

l*4500E+00 

2*0000E-0I 

MIX 


7*9500E-01 

0*0OOOE“Ol 

RECYCLE 

1*4600ET02 

l.OOOOE+OO 

MIX 


RECYCLE 2' 

0*00% OF * 3 






OUTFLOW TO 

♦ 2 INFLOW 

AIR TEMP 

= 2*1000EF02 F 

COMBUSTOR PRESSURE 

2*2100E+00 

ATM 

FUEL TEMP 

- 8*0000E+01 F 

LOWER HEATING VALUE 

l*8500E+04 

BTU/LBM 

SELECT AN 

OPTION BY number: 






•0- RUN WITH THIS DATA SET 


CHANGE 

NOMINAL COMBUSTOR 

PRESSURE 

-1- CHANGE 

AIR TEMPERATURE 

-5- 

CHANGE 

RECYCLE ELEMENT STATUS 

-2- CHANGE 

FUEL TEMPERATURE 

-6- 

CHANGE 

FLOW 

ELEMENTS STATUS 

-3- CHANGE 

LOWER HEATING VALUE -7- 

INSPECT 

SCHEMATIC MODEL 

LAYOUT 


OPTION? (0-7) 7 


MARK-II MODEL SCHEMATIC LAYOUT 


A F 
1 1 


A 

1 



A 

1 

^ ^ 






* 1 * 


* 2 * 

^ Z ^ 


* 4 

5KPSR*- 

— > 

•*MIX>f^-— 

»KPSR)1^- 


^^*:MIX* 


1 



1 

^ ^ 




PRESS <RETURN> TO CONTINUE 
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*** INPUT DATA 

FLOW AREA LENGTH FLOW INLET AIR INLET FUEL 

ELEMENT* <SQ*IN) (INCHES) TYPE (LBM/S) (LBM/S) 

1 l*4600E+02 4*0000E“01 PSR l*3450E+00 9>7200E-02 

2 l*4600E+02 1,0000E-01 MIX 6^9A00E-01 O.OOOOE-Ol 

3 1 ♦4600E+02 .1 *5000E+00 PSR 0*0000E~01 0,0000E~01 

4 l*4500E+00 2*0000E-01 MIX 7.9500E-01 O.OOOOE~01 

AIR TEMP = 2.1000E+02 F COMBUSTOR PRESSURE == 2*2100E+00 ATM 

FUEL TEMP 8*OOOOE+Ol F LOWER HEATING VALUE - l*8500E+04 BTU/LBM 

SELECT AN OPTION BY NUMBER: 

-0- RUN WITH THIS DATA SET -4- CHANGE NOMINAL COMBUSTOR PRESSURE 

-1- CHANGE AIR TEMPERATURE CHANGE RECYCLE ELEMENT STATUS 

-2~ CHANGE FUEL TEMPERATURE -6- CHANGE FLOW ELEMENTS STATUS 

-3- CHANGE LOWER HEATING VALUE -7- INSPECT SCHEMATIC MODEL LAYOUT 

OPTION? (0-7) 6 


TO CHANGE DATA, ENTER THE CODE AND ELEMENT (IE -B2* FOR *2'S LENGTH) 

CODE -A- CODE 'B- CODE “C* CODE Tr- 

FLOW AREA LENGTH FLOW INLET AIR 

ELEMENT* (SQ*IN) (INCHES) TYPE (LBM/S) 

1 1.4600E+02 4*0000E-01 PSR 1 *3450E+00 

2 l*4600E+02 1,0000E“01 MIX 6.9600E-01 

3 1 ♦ 4600E+ 02 1 ♦ 5000ETOO PSR 0 * OOOOE-01 

4 1»4500ET02 2»0000E-01 MIX 7*9500E-01 

CODE "F- “ REMOVE ENTIRE FLOW ELEMENT 

CODE 'G* - ADD ANOTHER FLOW ELEMENT (NO NUMBER NEEDED) 

(BLANK ENTRY RETURNS TO MAIN MENU) 

ENTER CODE, ELEMENT* A4 

ENTER CROSS-SECTIONAL AREA OF ELEMENT ♦ 4 (SQUARE INCHES): 2»00E2 
ENTER CODE, ELEMENT* A4 

ENTER CROSS-SECTIONAL AREA OF ELEMENT * 4 (SQUARE INCHES): 1,45E2 
ENTER CODE, ELEMENT* 


INPUT DATA 


Fl-OW 

AREA 


LENGTH 

FLOW 

IN’LET AIR 

INLET FUE 

ELEMENT* 

(SQ* IN) 


(INCHES) 

TYPE 

(LBM/S) 

(LBM/S) 

1 

1»4600E+02 


4,0000E 01 

PSR 

1 ,34S0E+00 

9.7200E-02 

r> 

l*4600E+02 


1*0000E-01 

MIX 

6*9600£-01 

O.OOOOE-01 

3 

1*4600E+02 


1 tSOOOEiOO 

PSR 

0*0000E-01 

O.OOOOE-01 

4 

l*4500E+00 


2.0000E--01 

MIX 

7.9500E-01 

0,OOOOE-Ol 

RECYCLE 

l,4A00E+02 


1*0000ET00 

MIX 

RECYCLE 20 
OUTFLOW TO 

.00% OF * 3 
* 2 INFLOW 

AIR TEMP = 

2*1000ET02 

F 

COMBUSTOR 

PRESSURE 

2*2100E+00 

ATM 

FUEL TEMP 

8.0000E+01 

F 

LOWER HEATING VALUE == 

1*8500ET04 

BTU/LBM 


SELECT AN OPTION BY NUMBER: 

-0- RUN WITH THIS DATA SET 
-1- CHANGE AIR TEMPERATURE 
-2- CHANCE FUEL TEMPERATURE 
-3- CHANGE LOWER HEATING VALUE 
OPTION? (0-7) 2 

ENTER FUEL TEMPERATURE (DEG FAHRENHEIT): 90* 


-4- CHANGE NOMINAL COMBUSTOR PRESSURE 
-5- CHANGE RECYCLE ELEMENT STATUS 
-6- CHANGE FLOW ELEMENTS STATUS 
-7- INSPECT SCHEMATIC MODEL LAYOUT 


CODE -E* 
INLET FUEL 
(LBM/S) 
9,7200E-02 
0,OOOOE-Ol 
0*OOOOE-Ol 
0»OOOOE-O:l 


193 



JK** INPUT DATA 


FLOW 

AREA 

LENGTH 

FL.OW 


INLET AIR 

INLET FUEL 

ELEMENT* 

(SQ*IN) 

(INCHES) 

TYPE 


(LBM/S) 

(LBM/S) 

1 

l,4600E+02 

4,0000E--01 

PGR 


l*34fJ0E+00 

9*7200E-02 

2 

1.4600E+02 

1*0000E“01 

MIX 


6.9600E-01 

O.OOOOE'OL 

3 

1.4600E+02 

l*5000E+00 

PBR 


0*0000E“01 

0*0000E-Ol 

4 

1*4500ET00 

2*0000E-01 

MIX 


7*9500E“01 

O.OOOOE-01 

r^ECYCLE 

l*4600E+02 

l^OOOOETOO 

MIX 


RECYCLE 20.00% OF ♦ 3 






OUTFLOW TO 

♦ 2 INFLOW 

(MR TEMP =: 

2.1000E+02 F 

COMBUSTOR PRESSURE 

2*2100E+00 

ATM 

FUEL TEMP - 

9»0000E+01 F 

LOWER HEATING VALUE = 

1 .8500E+04 

BTU/LBM 

SELECT AN OPTION BY NUMBER: 






-0- RUN WITH THIS DATA SET 

.-.4- 

CHANGE 

NOMINAL COMBUSTOR 

PRESSURE 

-1- CHANGE 

AIR TEMPERATURE 

-5- 

CHANGE 

RECYC 

LE ELEMENT STATUS 

-2- CHANGE 

FUEL TEMPERATURE 

-6- 

CHANGE 

FLOW 

ELEMENTS STATUS 

-3- CHANGE 

LOWER HEATING VALUE -7- 

INSPECT 

SCHEMATIC MODEL 

LAYOUT 


OPTION? (0-7) 0 

RUN WITH THIS DATA SET ARE YOU SURE? (Y/N) Y 


t** output: *** 

TEMP AND COMPOSITION FOR MARK-II MODEL: 

PRESSURE = 2,21 ATM OVERALL E.R. = 0.4995 

FUEL FLOW = 349,92 LBM/HR AIR FLOW = 2.84 LBM/SEC 

20.0 PCT. OF FLOW RECYCLED FROM OUTLET OF 
ELEMENT NO. 3 TO INLET OF ELEMENT NO, 2 

OTHER output: SELECT BY NUMBER - 

0- OUIT /BEGIN NEXT INPUT DATA SET 

1- INSPECT MOLE FRACTIONS 

2- INSPECT MASS FRACTIONS 

3- INSPECT EMISSIONS INDEX 

4- INSPECT MISC. OTHER OUTPUT 


5- INSPECT 
OPTION? (O' 

CURRENT INPUT 
-4) 1 

BATA SET 





MOLE FRACTIONS X(I) 






SPECES 


unt 

f2» 

5»C3* 

)»C43|£ 

EQL 

RCY 

C12H23 


7*484D-05 

4*070D-05 

2*1540-07 

1*5720-07 

2*8860-19 

2*1540- 07 

C2H2 


l*080D-02 

6*32211-03 

2*3000-03 

1*6790-03 

2*8860-19 

2*3000-03 

C2H3 


9*7a2D-04 

5*3270-04 

6*5580-06 

4*7850-06 

2*8860-19 

6*5580-06 

C2H4 


7*447D“04 

4*054D“04 

4*6250-06 

3*3750-06 

2*8860-19 

4*6250-06 

CH20 


5»914D“05 

3*2980-05 

4*3150-06 

3*1490-06 

1*6960-18 

4*3150-06 

CH3 


5*570D“05 

3*0920-05 

3*3320-06 

2*4310-06 

2*8860-19 

3*3320-06 

CH4 


3*81511-12 

2*5900-12 

2*6930-12 

1*9650-12 

2*8860-19 

2*6930-12 

CO 


1 *428D-02 

8*3350-03 

2*9020-03 

2*1180-03 

1*6940-06 

2*9020-03 

C02 


9*676Ei“02 

6*9700-02 

8*6160-02 

6*2880-02 

6*6580-02 

8*6160-02 

HCO 


l*064H-04 

5*8480-05 

3*3660-06 

2*4560-06 

5*7520-15 

3*3660-06 

H 


l*246D-03 

7*2750-04 

2*5540-04 

1*8640-04 

1*8140-08 

2*5540-04 

H2 


2*892D-03 

1*6810-03 

5*5690-04 

4*0640 04 

5*7610-07 

5*5690-04 

H20 


1*108D-01 

7*7160-02 

8*5410-02 

6*2330-02 

6*5690-02 

8*5410-02 

H02 


2*86111-04 

l*724ti-04 

8*5720-05 

6*2550-05 

7*4800-08 

8.5720-05 

N 


3.167D-07 

1.7510-07 

3*4790-08 

2*5390-08 

1*6020-13 

3*4790-08 

NO 


3*253D-05 

4*9410-05 

5*8140-05 

4*2430-05 

1*2010-03 

5*8140-05 

N02 


l*374D-07 

2*8790-07 

4*4420-07 

3*2410-07 

5*2850-06 

4*4420-07 

N2 


7*249D-01 

7*4690-01 

7*5120-01 

7*6170-01 

7*6340-01 

7*5120-01 

0 


1*53111-03 

1*0220-03 

9*5940-04 

7*0010-04 

1*9120-06 

9*5940-04 

OH 


4*489D-03 

2*9650-03 

2*6510-03 

1*9350-03 

7*1520-05 

2*6510-03 

02 


3*000D-02 

8*3860-02 

6*7440-02 

1*0600-01 

1*0100-01 

6*7440-02 

PRESS <return: 

> TO CONTINUE 
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output: 


TEMP AND COMPOSITION FOR MARK-II MODEL: 

PRESSURE ^ 2f21 ATM OVERALL E.R* = 0,4995 

FUEL FLOW = 349,92 LBM/HR AIR FLOW = 2,84 LBM/SEC 

20,0 PCT, OF FLOW RECYCLED FROM OUTLET OF 
ELEMENT NO, 3 TO INLET OF ELEMENT NO. 2 

OTHER output: SELECT BY NUMBER - 

0- QUIT / BEGIN NEXT INPUT DATA SET 
1“ INSPECT MOLE FRACTIONS 

2- INSPECT MASS FRACTIONS 

3- INSPECT EMISSIONS INDEX 

4- INSPECT MISC, OTHER OUTPUT 

5- INSPECT CURRENT INPUT DATA SET 
OPTION? (0-4) 4 


ELEMENT NO, *1* *2# *3* 

ELEMENT TYPE PSR MIX PSR 

EQUIV RATIO l,053D+00 A,941D-01 6,941D-01 
RES TIMEf SEC 5.647D-04 9,558D-05 1.261D-03 
AREA, SQ IN l,460D+02 1.460D+02 l,460D+02 
VELO, FT/SEC 5.903D+01 8,719D+01 9.909D+01 
FLOW, LBM/SEC l,442D+00 2,673D+00 2,673D+00 
AXIAL LOC, IN 4,000D-01 5,000D-01 2.000D+00 
ENTH , BTU/LBM-2 ♦ 219D+01-4 , 426D+00-4 , 427D+00 
TEMP EFF 8.192-01 7.850D-01 9.317D-01 

TEMP, DEG F 3.105D+03 2,400D+03 2.810D+03 
PRESS <R£TURN> TO CONTINUE 


*4* 

MIX 

4.995D-01 
1.879D-06 
1 ,45011+00 
8.872D+03 
2.933D+00 
2.200D+00 
5.530D+00 
9.260D-01 
2,192D+03 


EQL 

EQL 

4.995D-01 

O.OOOD-01 

0,000D-01 

O.OOOD-01 

2.933D+00 

2.200D+00 

5,530D+00* 

1 .OOOD+00 

2.350D+03 


output: 

TEMP AND COMPOSITION FOR MARK-II MODEL: 

PRESSURE == 2,21 ATM OVERALL E.R. = 0.4995 

FUEL FLOW == 349,92 LBM/HR AIR FLOW = 2.84 LBM/SEC 

20.0 PCT. OF FLOW RECYCLED FROM OUTLET OF 
ELEMENT NO. 3 TO INLET OF ELEMENT NO. 2 

OTHER output: SELECT BY NUMBER - 

0- QUIT / BEGIN NEXT INPUT DATA SET 

1- INSPECT MOLE FRACTIONS 

2- INSPECT MASS FRACTIONS 

3- INSPECT EMISSIONS INDEX 

4- INSPECT MISC. OTHER OUTPUT 

5- INSPECT CURRENT INPUT DATA SET 
OPTION? <0-4) 3 


RCY 
RCY-MIX 
A.941D-01 
4.205D-03 
1.460D+02 
1 .982D+01 
5*346D-01 
4,304D+10 
4,380D+00 
9.317D-01 
2,81011+03 
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EMISSION 

INDEX EI(I) , GM 

i I/KG FUEL 




SPECES ■! 

■ *1* 


)|c3* 


EQL 

RCY 

C12H23 

6.545D+00 

6.552D+00 

3*4490-02 

2*7590-02 

5*0490-14 

6*8970-03 

C2H2 

l*470D+02 

l*5B4D+02 

5*7310+01 

4*5850+01 

7*8580-15 

1*1460+01 

C2H3 

1*38311+01 

1*386D+01 

1*6970-01 

1*3580-01 

8*1620-15 

3*3940-02 

C2H4 

1*092D+01 

1*095D+01 

1. 2410 01 

9*9300-02 

8*4660-15 

2*4830-02 

CH20 

9*282D-01 

9*530D-01 

1*2400-01 

9*9130-02 

5*3250-14 

2*4790-02 

CH3 

4*377D-01 

4*47311-01 

4.7930-02 

3*8340-02 

4*5370-15 

9*5860-03 

CH4 

3. 199D-08 

3*999D-08 

4.1340-08 

3.3070-08 

4*8410-15 

8*2670-09 

CO 

2*091D+02 

2*246D+02 

7*7780+01 

6*2220+01 

4*9630-02 

1*5560+01 

C02 

2*226D+03 

2*952D+03 

3.628Dt-03 

2*9020+03 

3*1560+03 

7*2560+02 

HCO 

l*A14D+00 

l*633D+00 

9 *3440-02 

7*4750-02 

1*7450-10 

1*8690-02 

H 

6*565D-01 

7*056D-01 

2*4630-01 

1*9710-01 

1*9120-05 

4*9260-02 

H2 

3*047D+00 

3*262D+00 

1*0740+00 

8*5920-01 

1.2140-03 

2*1480-01 

H20 

l*043D+03 

l*338D+03 

1*4720+03 

1*1780+03 

1*2380+03 

2*9440+02 

H02 

4*936D+00 

5*476D+00 

2*7070+00 

2.1660+00 

2*5820-03 

5*4140-01 

N 

2*31BD-03 

2*361D-03 

4.662D-04 

3*7300-04 

2*3460-09 

9*3250-05 

NO 

5*102D-01 

1*42711+00 

1*6690+00 

1*3350+00 

3*7700+01 

3*3380-01 

N02 

3*303D-03 

l*275D-02 

1*9550-02 

1*5640-02 

2*5420-01 

3*9100-03 

N2 

1.061D+04 

2*013D+04 

2*0130+04 

2*2380+04 

2*2360+04 

4.0270+03 

0 

1.281D+01 

1*574D+01 

1. 4690+01 

1*1750+01 

3.1980-02 

2*9370+00 

OH 

3*990D+01 

4*853D+01 

4*3140+01 

3.4520+01 

1*2720+00 

8*6290+00 

02 

5*018D+02 

2*5820+03 

2*0650+03 

3*5570+03 

3*3800+03 

4.1300+02 

PRESS <RETURN> TO CONTINUE 






output: 

TEMP AND COMPOSITION FOR MARK-II MODEL: 

PRESSURE = 2*21 ATM OVERALL E*R* = 0*4995 

FUEL FLOW = 349*92 LBM/HR AIR FLOW 2*84 LBM/SEC 


20*0 PCT* OF FLOW RECYCLED FROM OUTLET OF 
ELEMENT NO* 3 TO INLET OF ELEMENT NO* 2 


OTHER output: SELECT BY NUMBER - 

0- QUIT / BEGIN NEXT INPUT DATA SET 

1- INSPECT MOLE FRACTIONS 

2- INSPECT MASS FRACTIONS 
3“ INSPECT EMISSIONS INDEX 

4- INSPECT MISC. OTHER OUTPUT 

5- INSPECT CURRENT INPUT DATA SET 
OPTION? (0-4) 0 

no YOU WANT A HARD COPY OF THE INPUT AND OUTPUT DATA FOR THIS RUN? (Y/N) Y 


DO YOU WANT ANOTHER RUN? (Y/N) N 
QUIT ARE YOU SURE? (Y/N) Y 


— FORTRAN STOP 
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PRESENT STATUS OF PROJECT - OCTOBER 1982 


Aval lable now: 

IMPROVED CREK (-0D) 

Including CREK preprocessor routines, 
sample data dekc, 1976 CREK documentation 

EQLBRM and MARK21 

interactive programs 

CREK- ID 

with CREK preprocessor routines and 1982 paper 
FOR PRELIMINARY EVALUATION AND CRITICISM ONLY! 


ALL SOURCE CODES; ANS I -standard 1966 FORTRAN-IV, ASCII 
MEDIA: Listing, or 8” floppy disk, single density, 

IBM 37^0 formatted; DEC RXOl in RT-11 or FILES-11 
other media negotiable 
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FLAME RADIATION MEASUREMENTS 


R. W. Claus, F. M. Humenik and G. M. Neely 
NASA Lewis Research Center 
Cleveland, Ohio 44135 

BACKGROUND ; 

One of the fundamental combustion experiments being 
conducted here at the Lewis Research Center consists of a 
series of flame radiation measurements on a single JT 8 D tubular 
combustor assembled in a cylindrical housing having optical 
viewing ports at three axial stations. 

Preliminary flame radiation data were obtained with the 
JT 8 D tubular combustor during August 1979. Technical paper, 
TP-1722, which was published in February 1981 presents some of 
the spectral radiation data which shows the relative radiance 
at three axial locations in the combustor for a few test 
conditions generally representative of low power settings. 
Supplementary radiant heat flux data, which were measured with 
thermopile type heat flux transducers at two axial stations 
were presented at the LeRC Symposium "Aircraft Research and 
Technology for Future Fuels", April 1980 (NASA Conference 
Publication 2146, pp. 153-160). All of these data were 
obtained with Jet A and ERBS fuels, and were limited to 
relatively low operating pressures and relatively mild 
operating conditions. Thus, a continuation and expansion of 
the previous experimental arrangement is required to include a 
broader range of operating conditions (P 3 , T3, f/a, fuel 
type) so that flame radiation characteristics can be mapped 
more fully. 

Since combustion gases, predominantly carbon dioxide and 
water vapor, emit and absorb radiant energy in discrete 
spectral bands, calculations of flame temperature and soot 
concentration can be performed using the experimentally 
measured spectral radiance profiles. In addition, these 
spectral radiation data can be used for possible correlations 
with local liner temperatures and/or the radiant heat flux 
measurements. Characteristic differences attributed to fuel 
property variations may also be evident from the spectral 
radiance signatures. Thus, there is considerable potential for 
obtaining very valuable data which could improve our 
understanding of the effects of operating conditions, fuel 
type, and other variables on total flame radiation. 
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RESEARCH OBJECTIVES: 


1. Improve and develop flame emissivity correlations. 

2. Evaluate effects of combustor pressure, inlet air 
temperature, and reduced fuel hydrogen content on 
flame radiance. 

3. Correlate increases in flame radiance with increases 
in average liner temperature. 

SIGNIFICANCE - APPLICATIONS : 

1. Establish a comprehensive flame radiation model. 

2 . Calculate flame temperature and soot concentration. 


PRELIMINARY RESULTS: 

Spectral and total flame radiation measurements have 
exhibited the following characteristics: 

1. Radiant heat flux increases substantially with rising 
combustor inlet air pressure. 

2. Fuel atomization characteristics can have substantial 
effects on radiant heat flux. 

3. Tests with ERBS fuel indicate that a reduction in fuel 
hydrogen content produces a significant increase in radiant 
heat flux primarily at low combustor pressures. At high 
combustor pressures (14-20 atm.) fuel effects are less 
distinct. 


REFERENCES: 


1. Claus, Russell W. : Spectral Flame Radiance from a 

Tubular-Can Combustor. NASA TP-1722, 1981. 

2. Aircraft Research and Technology for Future Fuels, NASA 
Conference Publication 2146, pp . 153-160; April 1980. 
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FLAME RADIATION MEASUREMENTS 


OBJECTIVE: IMPROVE LINER THERMAL ANALYSIS 


BENEFITS: t REDUCE DEVELOPMENT COST 

• IMPROVE COOLING AIRFLOW DISTRIBUTION 

• INCREASE LONG TERM DURABILITY OF LINER 


APPROACH; EXPERIMENTALLY ESTABLISH THERMAL MODEL PARAMETERS 

THROUGH DIRECT MEASUREMENT OF: 

1. SPECTRAL FLAME RADIANCE 

2. TOTAL RADIATIVE HEAT FLUX 



r THERMOCOUPLE 
■ RAKES I8S 


heat FLUX 

TRANSDUCERS 


IGNITOR 


r CAS SAMPLE 
' PROBES (41 


25 ATM 

870 K 
10 kg/sec 


COMBUSTOR 


LINER 

THERMOCOUPLES 
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CO/VIBUSTOR 


HWFLUX ^ 
TRANSDUCER 



EFFECT OF COMBUSTOR PRESSURE 


O 3 atm. Fuel nozzle #2 

□ 7 Jet A fuel 

A 14 

n 2 1 



.012 .016 .020 .008 .012 .016 .020 .008 .012 .016 


FUEL-AIR RATIO 


EFFECT OF FUEL ATOMIZATION 

O Fuel nozzle #1 Combustor pressure, 7 atm. 

O Fuel nozzle #2 Jet A fuel 



FUEL-AIR RATIO 



EFFECT OF FUEL TYPE 


□ JET A Combustor pressure. 7 atm. 

■ ERBS pygl nozzle *2 




FUEL-AIR RATIO 


EFFECT OF PRESSURE ON TOTAL FLAME RADIATION 
FOR TWO DIFFERENT FUELS. 



ZONE 3 

■% - 0.0 \(o 

T3 = £-3 3 K 

□ J"6r A 

M 07 : 7 .Lr: 2 




RADIANCE H/CM2 SR MICRON 


EFFECT OF PRESSURE ON SPECTRAL 
FLAME RADIANCE. 
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LEAN LIMIT PHENOMENA 


C. K. Law 

Northwestern University 
Evanston, Illinois 60201 

The influence of stretch and preferential diffusion on premixed flame 
extinction and stability have been investigated via two model flame 
configurations, namely the stagnation flame and the bunsen flame. The results 
are separately summarized in the following. 

(1) Extinction and Stability of Stretched Premixed Flames in the Stagnation 
Flow 

Using a counterflow burner and a stagnation flow burner with a 
water-cooled wall, the effect of downstream heat loss on the extinction 
of a stretched premixed flame has been systematically investigated for 
lean and rich propane/air and methane/air mixtures. Based on results of 
the concentration limits and flame separation distances at extinction, it 
is demonstrated that, in accordance with theoretical predictions, extinc- 
tion by stretch alone is possible only when the deficient reactant is the 
less mobile one. When it is the more mobile one, downstream heat loss or 
incomplete reaction is also needed to achieve extinction. A variety of 
unstable flame configurations have been observed; the mechanisms for 
their generation and sustenance are discussed. 

(2) Opening of Premixed Bunsen Flame Tips 

The local extinction of bunsen flame tips and edges of hydrocarbon/ 
air premixtures has been experimentally investigated using a variety of 
burners. Results show that, while for both rich propane/air and butane/ 
air mixtures tip opening occurs at a constant fuel equivalence ratio of 
1,44 and is therefore independent of the intensity, uniformity, and 
configuration of the approach flow, for rich methane/air flames burning 
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is intensified at the tip and therefore opening is not possible* These 
results substantiate the concept and dominance of the diffusional strati~ 
fication mechanism in causing extinction, and clarify the theoretical 
predictions on the possible opening of two-dimensional flame wedges. 

(3) Publications 

(a) **Lean“Lirait Extinction of Propane/Air Mixtures in the Stagnation- 
Point Flow,** by C, K, Law, S* Ishizuka, and M, Mizomoto, Eighteenth 
Symposium on Combustion , pp. 1791-1798 (1981) . 

(b) "Effects of Heat Loss, Preferential Diffusion, and Flame Stretch on 
Flame-Front Instability and Extinction of Propane/Air Mixtures," by 
S. Ishizuka, K, Miyasaka, and C. K. Law, Combustion and Flame , Vol. 
45, pp. 293-308 (1982). 

(c) "On the Opening of Premixed Bunsen Flame Tips,'* by C. K. Law, S. 
Ishizuka, and P. Cho, Combustion Science and Technology , Vol. 28, 
pp. 89-96 (1982). 

(d) "On Stability of Premixed Flame in Stagnation-Point Flow," by G. I. 
Sivashinsky, C. K. Law, and G. Joulin, Combustion Science and 
Technology , Vol. 28, pp. 155-159 (1982). 

(e) "An Experimental Study of Extinction and Stability of Stretched 
Premixed Flames," by S. Ishizuka and C. K. law, to appear in 
Nineteenth Symposium on Combustion , 1983. 
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OBJECTIVES 


• TO STUDY EFFECTS OF 

1. PREFERENTIAL DIFFUSION (Lef^l) 

2. AERODYNAMIC STRETCHING (FLOW NON-UNIFORMITY, 
UNSTEADINESS, AND FLAME CURVATURE) 

3. DOWNSTREAM HEAT LOSS 

• ON 

A. FLAME EXTINCTION 

B. FIAME-FRONT INSTABILITY 


METHODOLOGY 


1. PREFERENTIAL DIFFUSION EFFECTS STUDIED BY USING 

METHANE /AIR PROPANE /AIR 

LEAN l£ < 1 Le > 1 

RICH Le > 1 Le < 1 

2. AERODYNAMIC STRETCHING STUDIED BY USING STAGNATION 
FLOW WHICH HAS WELL-DEFINED VELOCITY GRADIENT 

3. DOWNSTREAM HEAT LOSS STUDIED BY USING 

(a) STAGNATION FLOW WITH WATER-COOLED SURFACE 

(b) SYMMETRICAL COUNTERFLOW 
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STAGNATION SURFACE 

/ / / // //////// /////////, 



Schematic of Stagnation-Point Flow Illustrating the 
Directions of Heat and mass Diffusion. 


EXTINCTION MECHANISMS 


Le > 1 Flames 

1. STRETCH ALONE CAN CAUSE EXTINCTION; DOWNSTREAM HEAT LOSS 
MINIMAL EFFECT 

2. INCREASING STRETCH EECREASES FIAME TEMPERATURE 

3. AT EXTINCTION, FLAME LOCATED AWAY FROM STAGNATION SURFACE 

4. AT EXTINCTION, DEFICIENT REACTANT COMPLETELY CONSUMED 


Le < 1 Flames 

1. INCREASING STRETCH INCREASES FLAME TEMPERATURE, THEREFORE 
STRETCH ALONE CANNOT CAUSE EXTINCTION 

2. EXTINCTION CAN BE ACHIEVED THROUGH 

(a) DOWNSTREAM HEAT LOSS, WITH FLAME AWAY FROM WALL 

(b) INCOMPLETE COMBUSTION WITH FLAME AT WALL 
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Various Flame Configurations for Propane/Air Mixtures 
in the Stagnation-Point Flow (See Publication No, b) 




Flame Separatedness at E; 
(a) Lean Methane/Air, (b 
(d) Rich Propane/Air (Se 


inct 

Rich 

Publ 











Ratio of the Extinction Concentration Limits with and without 
Downstream Heat Loss 


Author 

Method 


Definition 

% Methane 

% Propane 


Lean 

Rich 

Leon 

Rich 

Zabetakis 

Propagating flame 

(tube) 

Extinction 

5.0 

15.0 

2.1 

9.5 

Andrews and Bradley 

Propagating flame 

(vessel) 

Extinction 

^.5 

15.5 

- 

- 

Egerton and Thabet 

Flat flame 


Burning velocity 

5.1 

- 

2.01 

- 

Sorenson, Savage, 
and Strehlov 

Tent flame 


Cone angle 

^ 4.0 

15.0 

- 

- 

Yajnaoka and Tsujl 

Double flame 


Flame location 

I4.7 

15.3 

- 

- 

Ishizuka and Lav 

Binary flame 


Extinction 

U .8 

15.8 

2.0 

9.7 


Comparison of the Flammability Limits of Methane/Air and Propane/Air 
Mixtures Determined by Different Methods 
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[a] 


Schematic of Closed and Open Bunsen Flame Tips 



Set-Up of the Bunsen Flame Experiment 


216 



217 



(a) (b) 

Tip Intensification of Rich Methane/Air Bunsen Flame with 
Increasing Methane Concentration (See Publication No. c) 
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Tip Opening of Rich Propane/Air Bunsen Flame with Increasing 
Propane Concentration (See Publication No. c) 



VELOCITY CM/S) 


Fuel Concentrations at Tip Opening as Function of Flow Velocity for a Variety of Burners 
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COUPLING OF TRANSPORT AND CHEMICAL PROCESSES 
IN CATALYTIC COMBUSTION 


F.V. Bracco, C. BrunO/ B.S.H. Royce, D.A. Santavicca 
N. Sinha and Y* Stein 
Princeton University 
Princeton y New Jersey 08540 


Catalytic combustors have demonstrated the ability to operate effi- 
ciently over a much wider range of fuel air ratios than are imposed by the 
flammability limits of conventional combustors. Extensive commercial use 
however awaits further progress in the areas of i) the design of a catalyst 
with low ignition temperature and high temperature stability, ii) reducing 
fatigue due to thermal stresses during transient operation and, iii) the 
development of mathematical models that can be used as design optimization 
tools to isolate promising operating ranges for the numerous operating 
parameters. 

The catalytic combustion program at Princeton reflects these needs for 
further research. The current program of research involves the development 
of a two-dimensional transient catalytic combustion model and the develop- 
ment of a new catalyst with low temperature light-off and high temperature 
stability characteristics. 

A recently developed two-dimensional transient model has been 
used to study the ignition of CO/Air mixtures in a platinum coated cataly- 
tic honeycomb. This model includes coupling between the gas and substrate, 
radiative heat loss >to the outside and radial heat losses providing solu- 
tions for the two-dimensional temperature field in the substrate and the 
two-dimensional temperature, composition and velocity field in the gas. 

The model has been used to predict the transient and steady^ state using 
inlet conditions for which steady state experimental data is available. In 
the calculation the inlet conditions for velocity and temperature are set 
and the transient begins by injecting fuel. Results showing the transient 
behavior of the substrate and the gas are presented. This model is par- 
ticularly useful for predicting the temperature gradients in the substrate 
which in turn can be used to predict the thermal stresses in the catalyst 
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during the transient. Since transient measurements are not available the 
calculated steady state results are compared to steady state measurements. 
The predicted steady state catalyst temperature profile and the exhaust gas 
composition are found to agree very well with measurements over a range of 
inlet velocities and equivalence ratios. 

Operating temperatures of catalytic combustors often exceed 1500 K, 
resulting in unacceptably short lifetimes for standard catalysts. A new 
catalyst has been designed to exhibit long lifetime at high temperatures 
and adequate ignition characteristics at low temperatures. The catalyst is 
a modified perovskite based on La {Crg ^ 5 A£q^ 5 ) O 3 . This ceramic has been 
doped to make it electrically conductive and consequently it can be 
resistive ly heated to bring the catalyst up to the required light-off tem- 
perature. In addition^ platinum has been incorporated into the crystal 
structure to give improved low temperature light-off while having a low 
platinum vapor pressure at high operating temperatures. This catalyst can 
be used in powdered form, by washcoating it onto a high temperature ceramic 
substrate, or as sintered monolythic structures, e.g. plates. The advan- 
tage of using the catalyst in the form of plates is that they can be 
resistively heated to assist light-off. However, the catalyst is more 
readily available in powdered form and therefore the first tests with the 
new catalyst have been made with it washcoated on a honeycomb substrate. 

The substrate used was mullite, three inches long with 1/16 inch square 
cells. Results of experiments using a pure platinum washcoat and a 
perovskite powder with nominal 1 . 0 % (by weight) platinum washcoat are 
presented. The tests consist of establishing the light-off temperature and 
low temperature performance for several equivalence ratios and inlet velo- 
cities, after which the catalyst is "aged" for several hours at 1500 and 
then the light-off and low temperature experiments repeated. Preliminary 
results are very encouraging indicating very little change in surface acti- 
vity when used with propane fuel. Tests with hydrogen however indicate a 
poisoning effect which quickly reduces the activity of the catalyst. 

Results from these tests will be used as a guide in the selection of the 
optimum platinum content, in terms of both adequate low temperature light- 
off and minimal high temperature aging, to be used in the catalytic plates 
which will be produced and tested in the near future. 
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TRANSIENT CATALYTIC COMBUSTION MODEL 


• SOLVES COMPLETE COUPLED 2*D TRANSIENT CONSERVATION EQUATIONS IN 
BOTH SOLID AND GAS 

• PREDICTS TRANSIENT AND STEADY STATE, 2-D TEMPERATURE IN SUBSTRATE, 
2-D VELOCITY/COMPOSITION/TEMPERATURE IN GAS 

• INCLUDES RADIATIUU WlTrilN CATALYST BED AND FROM FRUNT/BALK TO 
UP/DOWN STREAM COMBUSTOR HALLS 

• INCLUDES RADIAL HEAT LOSSES (I.E- NON-ADIABATIC) 

• HAS BEEN USED TO PREDICT TRANSIENT AND STEADY STATE CO/AIR CATALYTIC 
COMBUSTION IN PLATINUM COATED HONEYCOMB CATALYSTS USING ONE STEP 
OVERALL REACTION FOR GAS AND HALL REACTIONS 


INLET CONDITIONS FOR TRANSIENT 

CALCULATION 

Velocity 

10 m/sec 

Temperature 

600 K 

Pressure 

110 KPa 

Fuel 

Carbon Monoxide 

Equivalence Ratio 

0.32 

Adiabatic Flame Temperature 

1680 K 

Water Mole Fraction 

5.^1 X 10-^ 

Mach Number 

0.02 

Reynold's Number 

290 
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LOCAL RADIATIVE LOSS/HEAT RELEASE 








* THE INITIAL HEAT RELEASE OCCURS NEAR THE ENTRANCE AT THE 
GAS-SOLID INTERFACE AND IS CONTROLLED BY HETEROGENEOUS 
REACTIONS. 

* LARGE SPATIAL AND TEMPERAL TEMPERATURE GRADIENTS OCCUR IN 
THE SOLID NEAR THE ENTRANCE CONTROLLED MOSTLY BY THE 
AVAILABILITY OF FUEL. 

* THE TEMPERATURE OF THE SOLID NEAR THE ENTRANCE ACHIEVES ALMOST 
ITS FINAL VALUE BEFORE SIGNIFICANT HEATING OF THE BACK. 

* HETEROGENEOUS REACTIONS AND THE GAS HEATED UP FRONT AND FLOWING 
DOWNSTREAM HEAT THE BACK OF THE SOLID. 

* THE OVERALL TRANSIENT TIME IS CONTROLLED BY THE THERMAL INERTIA 
OF THE SOLID AND BY FORCED CONVECTION. 

* RADIATION SIGNIFICANTLY INFLUENCES BOTH TRANSIENT AND STEADY 
STATE PARTICULARLY NEAR THE ENTRANCE. 

* THE OXIDATION OF CO OCCURS MOSTLY ON THE CATALYST AND SOON INTO 
THE TRANSIENT BECOMES DIFFUSION CONTROLLED. 
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I ) Two Objectives 

a) Low Temperature Light-Uff; Pt in Crystal Structure 

b) High Temperature Stability, Particularly Against Loss of Pt 


m) Perovskites, ABO3, Chosen as Class of Material to be Employed 



CRYSTAL STRUCIURE 
Approximately Cubic 
Easy Substitution of Impurities 
At a and B Sites 
To Control Properties 


1) Perovskites have been evaluated as MHU materials, and have proven high 
temperature endurance and catalytic activity 

2) The material is electrically conductive and may be heated by applied 
voltage, enabling precise light-off control 

3) Platinum dopant added for low temp, light-off capability; the Pt is 
chemically bound in the crystal structure for stability and endurance 

4) Honolythic design overcomes drawbacks of inert substrate + active 
coating catalysts 
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Ill) Compositions of Materials Fabricated 

a) (Lao.95 Sro.05HCrq. 499 A*-0.5 PfO-OOll 03 

Made by Transtech and suppied both as a p.owder for slip coating and 

AS MONOLYTHIC PLATES AND SPACERS FOR PARALLEL PLATE REACTOR 

b) La [Mgq.05 Crq.ii? Alo.I|7 Ptq.OII 03 

Made by Prof- Harlan Anderson, University of Missouri, and supplied 
as a powder suitable for slip coating 

iv) Basis For These Choices 

a) Electrical Conductivity 



Data from 
Meadowcroft, D-b- 
Brit- J- App- Phys 2 
(1969) 122S 


B) High Temperature Stability 


? uof 

E 


801 




ATM: P, 


Al.)^ 

= 10* otm 
^ s id*’ otm 


1740*C 


40F 


'M 


H--.-;- 


0,2 


0,4 0.6 

AL CONTENT. X 


0.6 


1,0 


Data from 

Prof- H-0- Anderson 


Weight loss 
function of 
gas vixture 


(Hg/cfli^/h) .of 

A1 content: Tenpetature • Atm: Flowing 

(0.001 atm Oi^.l atm Ht04)}2); Flow rate • 1 
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V) TEST CONFIGURATIONS 

a) Slip Coated onto a substrate* Al205 is a preferred substrate for 

HIGH TEMPERATURE OPERATION* FOR LIGHT-QFF TESTS> MULLITE MAY BE 
USED* 

DISADVANTAGES Phase diagram problems; no direct electrical heating 

FOR LIGHT-OFF CONTROL* 

ADVANTAGES Permits rapid test of material under combustion 

CONDITIONS 

b) Monolythic Structures - Powders may be sintered to produce monolythic 
materials in the form of plates* Densities between and 
theoretical are achievable and provide adequate plate strength* 

DISADVANTAGES Thermal shock cracking of monolyths (need to control 
density and composition to minimize this*) 


£ ^2r 

u 


to j Cr,.. ^U),»5 


8b 111 

O * 


< 

5 

cz: 


Temperature range 

25 to 1000®C 


0.2 0.4 0.6 0.8 

At CONTENT, X 


1.0 


ADVANTAGES Electrical heating possible for reaction initiation* 
Good experimental access, especially for laser probes* 
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Temp. (K) Temp. (K) 



AGING OF PURE PLATINUM CATALYST 
PROPANE FUEL, EQUIVALENCE RATIO = 0.42 
INLET TEMP. = 688 K. INLET VELOCITY = 9.9 m/s 



PURE PLATINUM CATALYST 

HYDROGEN FUEL. EQUIVALENCE RATIO = 0.1 

INLET TEMP. = 590 K, INLET VELOCITY = 10.2 m/s 


231 



Temp. (K) Temp. (K) 



PURE PLATINUM CATALYST 

PROPANE FUEL, EQUIVALENCE RATIO = 0.28 

INLET TEMP. = 589 K, INLET VELOCITY =10.1 m/s 



PEROVSKITE WITH NOMINAL 1 % BY WEIGHT 
PLATINUM CATALYST 

HYDROGEN FUEL. EQUIVALENCE RATIO = 0.10 
INLET TEMP. = 591 K, INLET VELOCITY =10.1 m/s 
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Temp. (K) 



PEROVSKITE WITH NOMINAL 1% BY WEIGHT 
PLATINUM CATALYST 

PROPANE FUEL, EQUIVALENCE RATIO = 0.30 
INLET TEMP. = 698 K, INLET VELOCITY = 6.0 m/s 


233 




SPONTANEOUS IGNITION CHARACTERISTICS 
OF HYDROCARBON FUEL-AIR MIXTURES 


Arthur H. Lefebvre and Greg W. Freeman 
Purdue University 
Lafayette, Indiana 47907 

Although the subject of spontaneous ignition of liquid fuels has 
received considerable attention in the past, the role of fuel evaporation 
in the overall spontaneous ignition process is still unclear. A main 
purpose of this research is to carry out measurements of ignition delay 
times, using fuels of current and anticipated future aeronautical interest, 
at test conditions that are representative of those encountered in modern 
gas turbine engines. Attention is focused on the fuel injection process, 
in particular the measurement and control of mean fuel drop size and 
fuel-air spatial distribution. The experiments are designed to provide 
accurate information on the role of fuel evaporation processes in 
determining the overall ignition delay time. The second objective is to 
examine in detail the theoretical aspects of spontaneous ignition in 
order to improve upon current knowledge and understanding of the basic 
processes involved, so that the results of the investigation can find 
general and widespread application. 

Experimental 

The first phase of the experimental program, which has just 
commenced, will utilize gaseous fuels only; namely propane and vaporized 
Jet A fuel. Its purpose is to determine accurately ignition delay times 
under conditions where a value of zero can be assigned to the fuel 
evaporation time. The tests will be conducted over a range of pressures 
from 1 to 10 atmospheres at unvitiated air inlet temperatures up to 900K. 

In some tests the oxygen content or nitrogen content of the main air 
stream will be varied in order to ascertain the effect of O 2 /N 2 ratio 
on chemical delay time. The information provided in these tests will 
be utilized in developing the theoretical model for spontaneous ignition. 

The test rig employed is shown schematically in Figure 1. The test 
section is constructed from several lengths of 30 mm dia. stainless steel 
tubing to form a maximum total length of 1 meter. Multipoint fuel 
injection is employed to ensure that the fuel -air mixing time is always 
negligibly small in relation to the ignition delay time. The test 
section is water-cooled to eliminate the possibility of flashback occurring 
along the boundary layer adjacent to the duct wall. It is proposed to use 
the same test procedure as that used successfully by Spadaccini and Te Velde 
(1). It consists of establishing prescribed conditions of pressure and 
fuel and air flow rates and gradually increasing the inlet air temperature 
until autoignition occurs at the exit of the test duct. The occurrence of 
autoignition will be detected by thermocouple probes located at the duct 
exit and by photodetectors located at several positions in the test rig. 

Theoretical 

Considerable progress has been made in simplifying and shortening 
the procedures involved in the calculation of fuel spray evaporation times. 
In a previous study by Rao and Lefebvre [2] the evaporation time was 
obtained as: 
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( 1 ) 



1 - 

1 

1-f 

1-n 

1 , 

2/3 


8 k. In (1+B) 

a 

1 + 0.25 Rep, 

, ^32 

V. y 


where 


B = mass transfer number 

0^2 = Sauter mean diameter, m 


f = fraction of fuel in vapor form 
k, = air conductivity, J/ms K 

a 

c = air specific heat, J/kg K 
^a 

u' = rms value of fluctuating velocity 

3 

p, = air density, kg/m 

a 

3 

p^ = fuel density, kg/m 

n = fraction of fuel initially in vapor form 

2 

Vg = air viscosity, m /s 


A drawback to the above equation is that it assumes a constant value of 
B for the entire evaporation period. However, B is constant only for 
steady-state evaporation. During the heat-up period the value of B rises 
continuously from an initial low value up to the steady-state value. 

Thus, in practice, the value of B for insertion into Eq. (1) should be 
lower than the steady-state value by an amount that depends on the 
ambient air pressure and temperature, and on the boiling temperature of 
the fuel. 

To overcome this deficiency the mass transfer number, B, is replaced 
by the evaporation constant, X, using the relationship 

8 (k/Cp) In (1 + B) 

^ f n\ 


Substituting for In (1 + B) from Eq. (2) into Eq. (1) gives 


236 




where Aeff is an "effective" value of x which takes into account both 
convective effects and the reduced evaporation rate during the heat-up 
period. Values of have been calculated for liquid hydrocarbons 
ranging in normal boiling temperature from 420 to 560 K, evaporating in 
air at temperatures and pressures up to 2000 K and 2000 kPa respectively. 
Some typical values of X are given in Figure 2 for a gas pressure of 
2000 kPa. 
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Figure 1. Schematic Diagram of Test Rig. 
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UNSTEADY FLOW EFFECTS IN COMBUSTION SYSTEMS 


Malladi V. S\i>»baiah 
California Institute of Technology 
Pasadena , California 91109 

Some of the important results obtained recently at Caltech in Prof. Frank 
Marble's research group are presented here. The following research personnel 
at Caltech contributed to the work presented here: 

Professor Frank E. Marble 
Dr. Malladi V. Subbaiah 
Dr. Perry Norton 
Dr. Ann Karagozian 

A wide variety of combustion problems, including combustion instabilities 
and turbulent diffusion flames, appear to involve the entrainment and 
deformation of laminar flames by large vortex structures in the flow field. 
First, we examine some details of this process of laminar flame distortion 
by considering the interactions of time-dependent diffusion flames with two- 
dimensional vortices (Figs. 1 & 2) . For large values of the circulation P/D, 
the augmentation of the fuel consumption due to the vortex is proportional to 
pp2/3Di/3^ When the effects of finite chemistry are included, the increase 
of fuel consumption rate is governed by a time scale which depends on the 
chemical reaction time, t^^ (Fig. 3). If the products of combustion occupy 
more volume than the original reactants, the spiral flame will appear as an 
unsteady volume dilatation for times on the order of chemical time. This acts 
as an acoustic source and the interaction of a vortex and diffusion flame 
results in the generation of a pressure pulse. The peak pressure is 
proportional to ^D^' ^//tch ^i^d occurs after a delay proportional to the 
chemical time, t^,^. The results provide the fundamental structure for the 
mechanism of instability proposed by Rogers and Marble (1956) (Figs. A-7) . 

In the second part of the presentation, some results on the modelling of the 
non-steady combustion in burners for aircraft gas turbines will be given. 

The general aim of the work is to develop a one-dimensional model applicable 
to the NASA-Lewis Non-Steady Combustion Rig. In the present discussion, we 
emphasize the results of the non-steady flame model, which constitutes an 
important module in the over all description of the system. 

In an earlier investigation, a detailed model for the non-steady response of . 
a stabilized flame in a two-dimensional duct was developed using a flame sheet 
description. The results showed active response by the flame region at certain 
well-defined frequencies and suggest a possible mechanism of low frequency 
instability in a combustion system. 

The present model for combustion processes utilizes a two-phase combustion 
model which treats the flame zone as an ensemble of pockets of unburned gas 
and combustion products. The steady state flame development is shown in 
Figures 8-11. This steady state solution is perturbed by an imposed acoustic 
wave approaching the flame region from either the downstream or upstream 
region. The spectra of the reflection and transmission coefficients are 
shown in Figures 12-15. 
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Figure 1 An initially flat flame has been wound into a 
spiral by a vortex; in this case the vortex is 
located on the flame sheet. Locally, a piece of 
the flame initially of length 6 \ has been 
elongated to length 61 . Note the similarity 
present in the vortex structure, the vortex 
grows with time as ^/TT. 
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Duct Wau 



Figure 4 An idealized representation of the experimen- 

tal setup of Rogers and Marble (1956). A con- 
dition of steady burning is shown here. 



Figure 5 The same configuration as the previous figure, 

only under screeching conditions. The 
unstable oscillations corresponded to the 
transverse mode of the combustion chamber, 
thus the acoustic oscillations are from top to 
bottom in this figure. The flame sheets are 
wound up by vortices alternately shed from 
the top and bottom of the flameholder. 
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Figure 6 To calculate the two-dimensional acoustic 
field caused by the spiral flame, one imagines 
the spiral flame occupying the z axis in a 
three-dimensional region, and an observer at 
0\ a distance x from the spiral flame. The 
pressure pulse at O' can be obtained by 
superimposing the pressure pulses from a 
line of three- dimensional sources distributed 
along the z axis. 



Figure 7 The pressure pulse seen at a distance x from 

the vortex, vrhen x is large enough to lie in 
the far field, is given by equations 7.8. Here 
the pressure rise. P — Pp, made dimensionless 
by (a — 1) is plotted as a 

function of (f x/c)/tg.f^. 
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o 

o 



Figure 8 Distribution of Combustion Products, Steady Flow 

^■/a, - 0.02, - 4.5. 

o 



Figure 9 . Velocity Distribution of Combustible Mixture, 
Steady Flow. 
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o 



Figure lO Velocity Distribution of Combustion Products, Steady Flow, 


o 



Figure 11 Pressure Distribution, Steady Flow, 
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1 Coefficient, T^, Wave From Downstream. 
6 



Figure 14 Pressure Wave Reflection Coefficient, R^, Wave From Upstream. 



Figure 15 Pressure Wave Transmission Coefficient, T^, Wave From Upstream, 
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NUMERICAL MODELING OF TURBULENT COMBUSTION 


A. F. Ghoneim, A. J. Chorin, and A. K. Oppenheim 
University of California, Berkeley 
Berkeley, California 94720 


Our work in numerical modeling is focused on the use of the 
random vortex method to treat turbulent flow fields associated with 
combustion while flame fronts are considered as interfaces betv/een 
reactants and products, propagating with the flow and at the same 
time advancing in the direction normal to themselves at a prescribed 
burning speed. The latter is associated with the generation of 
specific volume (the flame front acting, in effect, as the locus of 
volumetric sources) to account for the expansion of the flow field 
due to the exothermicity of the combustion process. The model was 
applied to the flow in a channel equipped with a rearward facing 
step. The results we obtained revealed the mechanism of the forma- 
tion of large scale turbulent structure in the wake of the step, 
while it showed the flame to stabilize on the outer edges of these 
eddies^^^ . 

In the course of this year we concentrated upon three topics: 

1. Study of the fundamental aspects of the modeling technique 

2. Application of the model to the formation of a turbulent jet 

3. Development of novel concepts bringing forth the aerodynamic 
properties of turbulent flames 

On the first subject, three efforts have been successfully pur- 
sued. First is the development of the random element method, pre- 
senting the application of the random walk technique to model dif- 
fusion of energy. The concept of vortex sheets and vortex blobs is 
extended for this purpose to heat transfer sheets or temperature jump 
elements, and blobs of internal energy, that carry thermal gradients 
from heat conducting walls to the interior of the field and redistribute 
them by diffusion. Table 1 presents the fundmental idea of random 
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walk modeling of diffusion, exploiting the analogy between the 
Green junction of the diffusion equation and the probability density 
junction of a Gaussian random variable. Figures 1 and 2 show a 
comparison between the exact solution and the numerical solution for 
the problem of heating an infinite solid by an isothermal wall at 
y = 0, in terms of the heat flux profiles and the temperature 
profiles. Figures 3 and 4 depict the temperature profiles for a 
finite solid bounded by two isothermal walls, and an isothermal and 
an adiabatic wall, respectively. The extention of the procedure to 
handle two-dimensional diffusion is presented in Table 2. A one- 
dimensional approximation close to the walls is employed to imple- 
ment isothermal boundary conditions, while energy elements are used 
in the interior to diffuse the energy there. Figures 5 and 6 display 
the temperature profiles along the diagonal of a corner and a square, 
respectively, as compared to the corresponding exact profiles. 

The second fundmental problem we have solved is concerned with the 
generation of yorticity by the interaction between the pressure field and 
density gradients — a mechanism which is of importance to both flames 
and buoyancy effects. The motion of the temperature jump elements, 
governed by. the algorithm of the random element method, generates ele- 
ments of vorticity in the interior of the field while additional 
elements are generated by the no-slip condition at the walls. Table 3 
describes how random walk is used in conjunction with the principle of 
time splitting to solve the system of equations that describes natural 
convection flow over a vertical infinite isothermal wall. Figure 7 

shows a comparison between our solution (thick lines) and a finite 

(3) 

difference solution for two values of the Prandtl number' . 

The third problem in this category is that of flame propagation. 

In our original model an interface advection and propagation algorithm, 
based on the assumption that the flame is a jump in density, was used to 
propagate the flame. In order to treat flame propagation as that of 
a reacting surface governed by a finite rate of reaction, the problem 
was recast in terms of a reaction-diffusion equation in temperature. 

The algorithm of the random element method was extended to solve this 
problem by adjusting the strength of the temperature jump elements as 
they move according to the rate of reaction as described in Table 4. 
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The algorithm was tested by solving problems by the use of finite 
difference and finite element methods, demonstrating its capability 
of calculating flame propagation with proper accuracy. Figures 8 and 
9 show a comparison between the exact solution, evaluated for 
f(T) = T(l-T), and the numerical computations when the integration 
of the reaction part of the equation is done using a first-order 
Euler scheme in the first case, and an exact integral in the second 
case, respectively^^\ 

On the second subject, the random vortex method was applied to 

the problem of the formation of planar, two-dimensional turbulent jet 

at high Reynolds numbers. Figure 10 describes schematically the 

elementary processes of the random vortex method and how they are 

implemented to solve the convection-diffusion equation. The results, 

expressed in terms of the development of the vorticity field, are 

presented in Fig. 11. They reveal the formation of large scale 

turbulent eddy structures on both sides of the jet with a potential 

core inside. Few jet widths downstream , the two layers start to 

interact and the flow becomes dominated by the pairing of eddies from 

both sides. The turbulent eddies grow by entraining the non-turbulent 

fluid, while their trajectories become more and more convoluted as a 

result of the interaction between positive and negative vorticity. 

These pictures display a remarkable resemblence to experimental photo- 

(51 

graphs obtained by Dimotakis, et al . ' ' across the plane of symmetry 
of an exi symmetric jet. 

The concept of flame aerodynamics has been developed on the basis 
of our numerical modeling studies, supported by experimental observa- 
tions of turbulent flames. The ultimate conclusion of these studies 
is the dominance of large scale eddies over the flow field, while the 
flame itself becomes established at the outer edges of these eddies 
and acts as a semi -permeable membrane encompassing the burnt gases. 

It is the interaction between. these eddies and the expansion associated 
with the exothermi ci ty of combustion process that produces the char- 
acteristic aerodynamic pattern that we are now studying^®^. 
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TABLE 1. FUNDAMENTAL IDEA 


Differential Equation 

3t 3^2 

Boundary Conditions 

4>(y.O) = 5(0) ; 4>(+ °°,t) = 0 

Constraint 

>|) = / 4> dy = 1 

— 00 

Formal Solution 

cj) = — ^ exp[-(y/a)^] ; a = v^2at 

/Zira 

Stochastic Solution 

E 6$. = 1 

y. (t + At) = y. (t) + 

E(n^ ) = 0 ; E(n^-^) = 2a At 

Local Sampling 

4, z {(y.y.) 

Global Sampling 

$ = E 5 <&. H(y-y.) 


Nomenclature 

H - Heavyside step function 
6 - Dirac delta function 

E - Expecte'd value 
n - Gaussian random variable 
a - Standard deviation 
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TABLE 2. A HYBRID SCHEME FOR TWO-DIMENSIONS 


DOMAIN 

BOUNDARY 
y < <5^ 

INTERIOR 

y > -Sg 

Differential Equation 

3T _ a^T 

“ 9y2 

at ' 

Boundary Condition 

T = 1 

T = 6(y-6^) 

Diffusing Element 

6T. 

56i 

Coupling 

6e^. = 6T^ * (y^'-Sg) 

6T, = ± 6e,/6^ 

Stochastic Solution 

r^. (t+At) = r^(t) + 

r^(t+At) = r^(t) + 

Sampling 

Z 6e. 
^ “ 6^-dx 

T - E 6e, 6(r-r,) 


^2 = 2a ; thickness of 1-D diffusion layer 
6A^. = area element 


254 










TABLE 3. PRESSURE-DENSITY INTERACTION 



n - number of temperature Jump elements 
N = Number of vortex elements 
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TABLE 4. REACTION-DIFFUSION EQUATION 


DIFFERENTIAL EQUATION 
INITIAL CONDITION 
BOUNDARY CONDITIONS 

_ 9 T ^ 

T(x,0) = Tq(x) 

T(+«^,0) = 0 ; T(-«>,0) = I 

SOLUTION 

FRACTIONAL STEPS 

REACTION 

DIFFUSION 

T(nAt) = [R(At) D(At)]" T^ 

w) 

dj _ dh 
3x^ 

DETERMINISTIC 

6T.(t + At) = 6T.(t) + f(T.) At 

STOCHASTIC 

T = 1 6T. H(x-x^) 


x^. (t + At) = x^ (t) + n^- 


E[n] = 0 ; E[n^] = 2 At 
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COMBUSTOR FLAME FLASHBACK 


James S* T'ien and Margaret Proctor 
Case Western Reserve University 
Cleveland, Ohio 44135 


Objective : Analytically model and conduct basic experimental tests to study 

the fundamentals of flame flashbacks in combustors. 

Background : Flashback has been a recurrent problem with the present Chrysler 

automotive gas turbine combustor. It can be a potential problem for 
other types and applications of advanced combustors as well. An improved 
understanding of the phenomenon would lead to improved design techniques 
to avoid its occurrence. 

Approach: It is proposed to model, test, and develop the fundamental 

conditions and process by which flashback occurs on combustor wall 
surfaces. An atmospheric rig consisting of a small dump combustor with a 
premixing channel would be required. Operation would be on a gaseous 
fuel. The primary variables would be inlet air and wall temperatures (up 
to 1500® F), boundary layer thickness, gas stream velocities, and 
controlled pressure disturbance level in the combustor. Both 
steady-state and transient tests will be conducted. 

Status : A steady-state flame flashback model is in existence. A transient 

model is being considered. The experimental rig is being designed. The 
experiment will be performed at NASA Lewis Research Center by a Case 
Western Reserve University student. 
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